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ABSTRACT 


A  study  was  conducted  to  determine  the  origin  of  features  appear¬ 
ing  in  a  Landsat  image  (1724-14472;  17  July  1974)  of  the  waters  off 
Cape  Cod  by  multitemporal  and  multispectral  analysis  techniques. 
Acquisition  of  ancillary  data  supplemented  imagery  from  five  different 
sensors,  including  visible,  near  infrared,  and  microwave,  representing 
fifteen  discrete  times  over  a  five-year  period.  Digital  data  for  the 
Landsat  image  was  evaluated  multispectrally .  Various  water  reflectance 
models  were  applied  to  the  digital  data  to  determine  whether  the  features 
related  to  reflected  light  from  the  ocean  bottom,  turbidity,  or  sur¬ 
face  reflectance.  Parameters  were  estimated  from  historical  data. 

Results  of  model  calculations  were  compared  with  the  observed  radiances 
at  several  locations  in  the  scene.  Multitemporal  analysis  affirmed 
that  the  following  detailed  information  can  be  determined  from  space 
imagery:  Rapid  rates  of  coastal  change,  possible  uncharted  shoals, 
and  theorized  surface  manifestations  of  internal  waves.  In  all  cases 
multisensor  imagery  increased  the  reliability  of  the  derived  infor¬ 


mation.  Multispectral  analysis  provided  additional  evidence  for 
attributing  variations  of  reflectance  in  deep  ocean  to  the  effects  of 
surface  reflection  and  also  appears  to  be  the  primary  cause  of  re¬ 
flectivity  of  patterns  over  the  shoal  areas  around  Nantucket  Island. 
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INTERPRETATION  OF  HYDROGRAPHIC  FEATURES 
IN  THE  WATERS  OFF  CAPE  COD 

by 

Charles  Leonard  Reed 

On  13  December  1976,  the  tanker  Argo  Merchant  was  stranded  on 
Fishing  Rip  in  the  Shoals  of  Nantucket.  Its  cargo  of  3  x  10^  liters 
of  crude  oil  spilled  into  the  sea,  focusing  American  attention  on  an 
environmental  disaster  of  unprecedented  magnitude.  Available  charts, 
aerial  photographs,  satellite  imagery,  and  environmental  background 
information  were  reviewed  to  help  assess  the  potential  damage.  Among 
the  available  data  was  Landsat  image  1724-14472  of  17  July  1974, 
showing  the  area  from  Glouchester,  Massachusetts,  to  Block  Island, 

Rhode  Island,  including  the  intervening  coastal  waters  and  the  Shoals 
of  Nantucket  (Figure  1)  (Short,  et  al.,  1976).  Recorded  in  this  image 
are  several  patterns  suggestive  of  bottom  features.  If  such  patterns 
can  be  linked  to  bottom  features,  Landsat  imagery  may  have  important 
applications  in  nautical  charting. 

OBJECTIVES  OF  THIS  STUDY 

The  study  described  in  this  report  was  conducted  to  determine  the 
origin  of  several  features  appearing  in  Landsat  image  1724-14472  of 
17  July  1974,  showing  the  Cape  Cod  area.  Specific  attention  was  given 
to  four  questions: 

1.  Are  the  features  detected  in  the  image  related  to  water  depth? 
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2.  Can  depch  information  be  derived  from  Landsat  data? 

3.  Do  the  features  detected  represent  hazards  to  navigation? 

4.  If  not  related  to  water  depth,  could  these  features  obscure 
or  mask  hazards  to  navigation? 

RATIONALE  FOR  THIS  STUDY 

The  compilation  of  charts  for  coastal  areas  and  shallow  waters  is 
perhaps  the  most  difficult  task  facing  the  hydrographic  charting 
community.  Conventional  precision  bathymetric  surveys  are  costly, 
both  in  terms  of  time  and  money.  The  coastline  and  bathymetry  are 
continually  changing  from  both  natural  forces  and  man-made  modifi¬ 
cations.  Available  charts  are  usually  inaccurate  due  to  dynamic 
changes  which  occur  during  the  time  it  takes  to  create  the  chart 
product.  All  marine  charts  issued  by  the  Defense  Mapping  Agency 
Hydrographic/Topographic  Center  (DMAHTC)  have  an  edition  date,  and 
should  be  continually  corrected  from  information  disseminated  through 
the  weekly  Notice  to  Mariners.  Such  chart  corrections  must  be  manu¬ 
ally  entered  on  the  charts  and  many  civilian  users  are  unaware  of  the 
numerous  changes  making  their  charts  obsolete,  especially  if  their 
chart  is  an  old  edition. 

Shipboard  bathymetric  coverage  of  most  ocean  areas  is  minimal  be¬ 
cause  ships  are  only  capable  of  linear  coverage,  requiring  a  grid  of 
transects  to  cover  the  survey  area.  Modern  multi- beam  sonar  systems 
are  of  little  help  in  shallow  waters,  for  they  do  not  provide  data 
for  much  more  of  the  ocean  bottom  than  do  conventional,  single  beam 
echo  sounders.  In  addition,  ships  conducting  hydrographic  surveys  in 
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uncharted  shallow  waters  are  subjected  to  greater  hazards  of  ground¬ 
ing  on  unexpected  shoals  or  reefs.  The  use  of  traditional  aerial 
surveys  is  hampered  by  adverse  weather  in  many  coastal  areas,  and  may 
be  additionally  limited  by  aircraft  range,  political  considerations, 
and  the  limited  availability  of  support  facilities.  If  shipboard 
survey  data  and  chart  compilation  resources  could  be  augmented  by  use 
of  Landsat  data,  better  charts  could  be  produced  faster,  and  at  less 
cost,  than  is  now  possible. 

While  the  ultimate  benefit  of  spacecraft  imagery  might  be  actual 
determination  of  ocean  depth,  especially  in  shallow  waters,  early 
detection  of  hazards  to  navigation  would  be  a  major  benefit.  Infor¬ 
mation  on  location  and  changes  in  underwater  obstacles  such  as  shoals 
could  be  used  to  alert  mariners  and  revise  existing  charts. 

Modern  nautical  chart  products  are  basically  concentrated  in  areas 
of  potential  military  interest,  and  along  major  shipping  lanes.  Be¬ 
cause  survey  resources  are  allocated  by  priority,  charts  of  many  areas 
of  the  world  are  so  out  of  date  or  non-existent  as  to  be  unresponsive 
to  the  changing  requirements  of  the  1980' s  (Dawson,  1980).  Most  charts 
of  the  world's  oceans  are  based  on  surveys  of  the  1800 's  (Hammack, 

1977).  These  charts  were  adequate  during  their  day  because  vessels  did 
not  exceed  6  meter  draft  and  sailed  upon  fairly  well  established  trade 
routes.  Today,  our  industrial  society's  appetite  for  raw  materials  has 
led  to  development  of  super  tankers  and  deep  draft  cargo  ships  to  carry 
increased  exports  from  developing  countries.  Modern  vessels  include 
the  largest  bulk  carriers  ever  devised  by  men,  with  some  over  360  meters 
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in  length  and  drawing  over  28  meters  (Cramer,  1975).  Former  marine 
charting  standards  are  not  adequate  to  meet  the  needs  of  these  new 
vessels,  and  the  number  and  cost  of  marine  accidents  is  increasing. 

In  an  economy  of  ever  rising  costs  and  inflation,  navigators  must  also 
seek  shorter  routes,  more  favorable  sailing  weather,  and  take  advan¬ 
tage  of  any  favorable  currents  to  minimize  fuel  cost  and  avoid  delays 
in  transit  (Dawson,  1980) . 

The  combined  elements  of  new  vessels,  new  ports,  and  new  trade 
routes  require  the  production  of  new  charts.  This  is  an  expensive 
undertaking  and  interest  in  satellite  technology  to  assist  in  this 
undertaking  is  high.  The  shortfall  in  chart  production  is  so  large 
that  new  technology  must  be  utilized.  To  meet  long-term  requirements 
for  precise  charting,  sensor  systems  such  as  Landsat  must  be  coupled 
with  limited  ship  surveys  to  find  a  feasible  solution  to  the  current 
charting  problem  at  a  reasonable  cost. 

CHARACTERISTICS  OF  THE  STUDY  SITE 

The  waters  around  Cape  Cod  represent  a  highly  dynamic  and  complex 
environment,  in  which  the  combined  forces  of  wind  and  water  act  upon 
the  surface  of  the  earth.  The  geologically  rapid  changes  in  land 
forms  due  to  erosion  and  aggradation  are  significant  to  topographic 
and  bathymetric  mapping,  and  represent  a  difficult  and  continuing  task 
of  map  revision.  The  cask  of  hydrographic  charting  which  includes 
both  shorelines  and  bathymetry,  is  more  difficult  than  mapping  of  the 


land  surface. 


The  waters  around  Cape  Cod  provide  a  textbook  example  of  multiple 


use  problems,  and  a  wide  variety  of  marine  charting  applications  are 
uniquely  represented  by  this  multiplicity  of  uses.  Oil  leases  have 
now  been  sold  on  the  Georges  Bank.  To  pump  oil  to  the  mainland  re¬ 
quires  routing  of  pipelines  and  cables.  Problems  of  oil  drilling  and 
pipeline  construction  are  compounded  because  the  area  is  one  of  the 
most  productive  fishing  areas  in  the  world.  Frequent  dredging  is  re¬ 
quired  to  allow  the  fishing  fleet  from  Chatham  to  sail  past  the  barrier 
island,  Nauset  Spit,  and  Monomoy  Island  to  fishing  grounds.  The  waters 
off  Cape  Cod  and  Nantucket  Shoals  represent  one  of  the  busiest  com¬ 
mercial  shipping  lanes  of  the  Atlantic  Ocean,  and  the  area  is  one  of 
the  major  maritime  graveyards  of  the  world.  More  than  145  shipwrecks 
have  been  recorded  in  the  area  since  World  War  II  (Green,  1979),  with 
perhaps  2,000  since  colonization  (Melham,  1975). 

Glacial  Origins 

The  unique  land  form  that  is  called  Cape  Cod  began  to  develop  at 
the  end  of  the  last  Ice  Age,  during  the  final,  or  Wisconsin  stage  of 
the  Pleistocene  Epoch  (from  about  50  to  70  thousand  years  ago  until  10 
thousand  years  ago).  As  the  glaciers  grew,  sea  level  declined  about 
90  meters  below  present  sea  level  and  the  area  we  now  know  as  Cape  Cod 
was  a  coastal  plain.  The  ice  spread  southward  with  three  separate 
lobes  converging  over  the  area  leaving  the  basic  deposits  which  formed 
Cape  Cod.  Deposition  occurred  in  two  stages.  At  their  furthest  ad¬ 
vance  southward,  the  three  lobes  of  the  ice  sheet  left  end  moraine 
deposits  which  formed  Martha's  Vineyard  and  Nantucket  Island  (Figure  2) 
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The  ice  then  retreated,  but  during  a  temporary  delay  in  this  retreat 
formed  the  south  shore  of  Cape  Cod  Bay.  Two  secondary  inter lobate 
moraines  were  formed  between  these  two  periods  of  glacial  equilibrium 
with  the  climate.  Both  were  somewhat  smaller  in  size  and  trend 
north- south  from  the  northernmost  point  of  Martha's  Vineyard  to  the 
west  side  of  Cape  Cod  Bay,  and  from  Nantucket  Island  along  the  east 
side  of  Cape  Cod.  No  remnant  of  the  latter  moraine  persists  due  to 
erosion  from  the  sea,  but  its  existence  can  be  inferred  from  the 
generally  smooth  appearance  of  the  sediment  from  the  ocean  to  Cape  Cod 
Bay.  The  "forearm"  of  Cape  Cod  was  an  outwash  plain  from  a  glacier  and 
moraine  deposit  on  the  Atlantic  side  (Strahler,  1966). 

Erosion  and  Accretion 

« 

After  the  major  glacial  deposits  were  laid  down  and  the  glaciers 
retreated,  sea  level  began  to  rise  to  the  present  level  and  a  second 
great  shaping  of  the  land  commenced.  Rough  edges  of  the  topography 
were  smoothed  by  wave  action  to  the  more  mature  shoreline  of  today. 

Erosion  and  accretion  are  very  rapid  in  unconsolidated  glacial 
sediments,  such  as  those  of  Cape  Cod.  The  high  energy  expended  on 
these  easily  eroded  deposits  by  wind,  tide,  and  sea  state  can  produce 
dramatic  changes.  Winds  are  normally  from  the  northwest  or  southwest, 
but  severe  storm  winds  are  often  from  the  northeast;  coloquially 
called  "nor'  easters."  Winds  from  the  northeast  encounter  no  ob¬ 
structions  across  the  entire  breadth  of  the  North  Atlantic,  and  large 
swell  waves  generated  from  the  large  fetch  of  the  ocean  batter  the 
shoreline  and  shoal  areas.  The  waves  strike  the  land  at  an  angle, 
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secting  up  littoral  currents  that  transport  sediment  along  the  shore, 
adding  sediment  to  offshore  spits,  and  sometimes  actually  closing  off 
estuaries.  Normal  tidal  range  is  between  3  and  4  meters  (Strahler, 
1966),  but  Cape  Cod  often  has  very  high  tides  associated  with  storms. 
Effects  of  tidal  surges  upon  the  land  form  are  pronounced  due  to  the 
wide  and  shallow  continental  shelf.  Intensely  rippled,  broad  tidal 
flats  characteristically  result  from  the  strength  of  the  tides. 

The  coastal  zone  is  a  dynamic  system  in  which  sediment  is  con¬ 
tinuously  in  motion  due  to  sea  state,  tides,  and  littoral  drift; 
therefore.  Cape  Cod's  shoreline  is  continuously  changing.  An  example 
of  long  term  coastal  processes  can  be  found  in  the  account  of  the 
first  recorded  shipwreck  in  the  New  World  at  Nauset  Beach,  Cape  Cod, 
Massachusetts,  in  1626.  The  Sparrow-Hawk,  a  12  meter  (40- feet) 
shallow  draft  vessel  set  sail  from  England  for  Virginia.  During  the 
crossing  of  the  Atlantic  the  captain  became  ill  with  scurvy,  and 
provisions  ran  low.  After  six  weeks  without  water,  "...  the  immigrants 
forced  the  crew  to  steer  between  southwest  and  northeast  in  hope  of 
finding  some  land,  what  so  ever,  and  caring  not"  (Anonymous,  1865; 
Onysko,  1978). 

The  ship  first  ran  aground  on  an  offshore  sand  bar  late  one  even¬ 
ing.  Rising  tides  the  following  morning  freed  it  to  sail  safely 
through  an  inlet  into  what  is  now  known  as  Pleasant  Bay.  A  violent 
storm  (perhaps  a  "nor1  easter")  subsequently  grounded  and  wrecked  the 
Sparrow-Hawk,  engulfing  it  in  sand.  After  the  storm,  the  inlet  mi¬ 
grated  to  the  south  and  coastal  processes  completely  buried  the  ship. 
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In  1863,  237  years  after  it  was  wrecked,  the  Sparrow-Hawk  was  uncovered 
by  yet  another  storm.  Since  1863,  the  inlet  has  again  migrated  away 
from  the  wreck  site. 

Short  term  affects  of  wind,  tide,  and  sea  state  upon  the  shoreline 
of  Nantucket  Island  have  been  described  by  Gutman,  et  al.,  (1979). 

The  island  can  lose  as  much  as  11  meters  to  the  sea  per  year,  depend¬ 
ing  upon  wind  and  wave  direction,  and  on  the  amount  of  energy  dissi¬ 
pated  from  storm  waves  by  the  sand  shoals  which  lie  to  the  north  and 
east  of  the  island. 

Sand  Shoals 

The  Nantucket  Shoals  represent  a  remarkable  depositional  topog¬ 
raphy.  The  ridges  and  trenches  of  the  shoals  trend  northwest- southeast, 
and  rise  to  within  4  meters  (2  fathoms)  of  the  surface.  The  Nantucket 
Shoals  deposition  has  been  compared  "...  to  the  Norfolk  Banks  of  the 
southern  bight  of  the  North  Sea,  or  the  Syrtis  Major  section  of  Mars" 
(Swift,  1978).  The  shoals  are  overprinted  by  a  second  pattern  of 
sand  "waves"  up  to  seven  meters  high  and  spaced  hundreds  of  meters 
apart.  These  features  are  clearly  visible  on  Sky lab- 3  S-190B  color 
photography  (Figure  3).  Similar  sand  ridges  on  Georges  Bank,  120  miles 
east  of  Cape  Cod,  were  studied  by  Stewart  and  Jordan  (1964)  after  the 
shoals  were  attributed  to  reshaping  of  glacial  deposits  by  strong  cur¬ 
rent  action  (Shepard,  et  al.,  1934).  Three  interesting  phenomena  were 
observed.  The  shoals  were  found  to  have  migrated  up  to  310  meters  over 
a  period  of  25  years;  the  tidal  pattern  was  found  to  move  rotationally , 
360°  in  12  hours,  rather  than  reciprocally;  and  a  change  in  level  of 


Che  water  surface  was  found  above  the  ridges  of  the  shoals.  The  rise 
in  water  level  was  attributed  to  an  overfall  created  when  tidal  cur¬ 
rents  move  normal  to  the  trend  of  the  ridges.  An  actual  rise  in 
the  sea  surface  can  be  detected  at  the  confluence  of  the  currents. 

Swift  (1978)  described  the  area  as  follows: 

"The  development  and  migration  of  the  shoals  are  a  response 
of  the  cohesionless  sand  substrate  to  the  intense  tidal  flow. 

The  ridge  channels  are  not  exact  replicas  of  each  other; 
some  shoal  north  and  some  shoal  south;  depending  on  whether 
they  service  primarily  the  south  moving  ebb  tidal  stream  or 
the  north  moving  flood  tidal  stream.  Sand  wave  crests  are 
not  truly  orthogonal  to  the  ridge  set;  sand  wave  crests 
tend  to  approach  ridge  crests  asymetrically ,  so  that  sand 
waves  are  convex  toward  the  direction  from  which  the  dominant 
tidal  stream  approaches,  and  the  residual  flow  moves  alter¬ 
nately  clockwise  and  counterclockwise  around  successive 
ridges . " 

McCave  (1971)  studied  similar  sand  shoaling  and  tidal  currents  in 
the  North  Sea  to  attempt  to  learn  more  about  the  depositional  mechanism 
of  sand  transported  by  tidal  currents.  Knowledge  of  the  growth  and 
migration  of  shoaling  could  be  very  useful  in  reducing  the  danger  to 
shipping  and  the  potential  environmental  dangers  inherent  in  transport 
of  oil  by  ships  or  pipelines  across  such  rich  fishing  grounds. 


Currents 

Major  currents  contribute  to  the  unique  marine  climate  of  Cape 
Cod.  The  Gulf  Stream,  which  carries  warm  water  from  the  tropics, 
turns  eastward  as  its  northward  flow  meets  the  cold  Labrador  Current 
flowing  southward  from  the  Arctic.  The  upwelling  of  nutrients  caused 
by  the  meeting  of  the  currents  promotes  the  growth  of  phytoplankton, 
the  primary  link  in  the  food  chain  for  the  fishing  grounds .  Approxi¬ 
mately  1.3  x  107  kilograms  of  commercial  fish  with  a  value  of  S 7 8 


million  are  taken  annually  from  waters  off  the  coast  of  Massachusetts 
(Fricke  and  Maiolo,  1978). 

Fog  and  marine  haze  are  characteristic  of  Cape  Cod  in  the  warmer 
months,  when  cold  waters  from  tidal  currents  are  warmed  over  shallows, 
such  as  the  ridges  of  the  shoals  (Ryder,  1979).  Fog  patches  up  to  10 
kilometers  across  are  common. 

Very  strong  tidal  currents,  with  velocities  up  to  three  knots, 
occurring  in  the  Cape  Cod  area  are  frequently  called  rip  tides.  Charts 
of  the  area  designate  many  sites  as  "rips,"  such  as  Fishing  Rip,  the 
site  of  the  Argo  Merchant  grounding.  The  along-shore,  or  littoral- 
drift,  current  resulting  from  the  oblique  angle  at  which  waves  strike 
the  shore,  is  a  dominant  force  in  the  transport  of  sediment  from  land 
erosion  and  subsequent  accretion. 

CHARACTERISTICS  OF  THE  DATA  SET 

This  study  was  stimulated  by  observations  of  potential  bottom 
features  in  Landsat  image  1724-14472  of  17  July  1974,  but  several 
additional  images  were  also  studied.  The  full  image  set  included 
records  from  14  dates,  representing  several  sensors  and  spectral 
bands  (Table  1) .  These  remote  sensor  data  were  supplemented  with 
data  on  weather  at  the  times  the  images  were  obtained  (Table  2), 
tidal  data  for  Nantucket  (Table  3),  and  tidal  currents  at  Pollock 
Rip  Channel  (Table  4). 

Most  of  the  remote  sensor  data  originated  with  a  sensor  on  one 
of  the  Landsat  series  of  satellites,  but  data  from  Sky  lab- 3  and 
Seasat-1  were  also  utilized.  In  each  case  the  record  represents  the 
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spatial  distribution  of  intensities  of  electro-magnetic  energy  in 
specific  spectral  bands  determined  by  sensor  configuration.  To  pro¬ 
vide  a  record  of  underwater  phenomena,  this  energy  must  travel  from 
its  source  --  the  sun  except  for  the  Seasat- 1  radar  data  --  to  the 
surface  of  the  water,  travel  downwards  through  the  water  column  to 
the  bottom,  be  reflected  off  of  the  bottom,  travel  upwards  from  the 
bottom  through  the  water  column  to  the  surface,  and  travel  from  the 
water  surface  to  the  sensor.  Thus,  each  record  represents  the  inte¬ 
grated  result  of  a  large  number  of  variable  factors  affecting  the 
"energy- flow- profile"  from  source  to  sensor  (Olson,  1962;  Colwell,  et 
al.,  1963).  Solar  energy  recorded  by  the  Landsat  and  Skylab  sensors 
included  two  components:  direct  solar  irradiation  and  indirect  irradi¬ 
ation  from  atmospheric  scattering  in  what  is  assumed  to  be  a  uniform 
Lambertian  sky  (Francis  and  Reed,  1979).  Scattering  of  energy  in  the 
atmosphere,  reflectance  from  the  sea  surface,  and  scattering  of  energy 
in  the  water  column  all  tend  to  increase  the  amount  of  energy  reach¬ 
ing  the  sensor  and,  therefore,  tend  to  reduce  the  contrast  between 
deep  water  and  bottom  signals  received  at  the  sensor.  The  final  re¬ 
sponse  of  the  sensor  is  a  product  of  the  sensitivity  of  the  sensor 
and  the  radiance  received: 


Vi  =  c(  X  )iL(  \  >i 


where , 


vi 

c(  X  ) 
L(  X  ) 


=  the  sensor  response 

=  the  spectral  sensitivity  of  the  sensor 
=  the  spectral  radiance  received 
=  the  spectral  band  of  operation  of  the  sensor. 
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The  spectral  radiance  received  at  the  sensor  includes  several 
components,  only  one  of  which  originates  from  bottom  reflectance. 
Specular  reflections  from  water  surfaces  are  corrsnon,  but  are  more 
limited  in  Landsat  data  than  in  conventional  photography  because  of 
the  narrow  acceptance  angle  (field-of-view)  of  the  Landsat  sensors. 
What  specular  reflections  do  occur  are  usually  the  result  of  sloping 
water  surfaces  of  waves  or  currents,  and  appear  as  "glitter"  in  image 
forms  of  the  data.  Within  the  water  column  of  the  ocean,  scattering 
is  minimal  and  little  of  the  scattered  energy  exits  the  water.  Thus, 
the  volume  scattering  of  the  water  column  is  usually  ignored.  The 
signal  reflected  from  the  bottom  is  strongly  affected  by  attenuation 
in  the  water  column.  This  attenuation  is  proportional  to  the  clarity 
and  type  of  water,  and  to  the  length  of  the  water  column.  Scattering 
in  the  atmosphere  adds  a  non- directional  radiance  to  the  total  signal 
reaching  the  sensor.  Thus,  the  total  radiance  recorded  by  a  sensor 
over  water  is  composed  of  three  primary  components,  and  can  be  ex¬ 
pressed  as  follows: 

Ltotal  =  Lp  +  Ls  +  Lb 

where,  Lp  =  the  atmospheric  path  radiance  due  to  scattering 
Ls  =  the  surface- ref lected  radiance 
L^  =  the  bottom- reflected  radiance 

both  Lp  and  Ls  are  assumed  to  be  homogeneous  for  a  scene.  Lj,  is  a 
result  of  bottom  reflectance  and  water  attenuation.  The  atmospheric 
path  radiance,  Lp,  and  surface- ref lected  radiance,  Ls,  contain  no  in¬ 
formation  regarding  water  depth,  and  are  therefore  considered  "back¬ 
ground"  radiance  which  must  be  subtracted  from  the  total  radiance  for 


-13- 


water  depth  calculations.  The  sum  of  these  two  factors  is  the  radiance 
of  deep  ocean  water,  L^.  Unfortunately,  there  is  nothing  in  the  signal 
itself  to  identify  what  part  of  the  signal  reaching  the  sensor  origi¬ 
nated  as  path  radiance,  as  surface- ref lected  radiance,  or  as  bottom 
reflectance.  Accuracy  of  bathymetry  extracted  from  Landsat  would  be 
degraded  by  erroneous  interpretation  of  signals  received  from  haze 
over  the  water,  foam  on  the  water,  or  vegetation  in  the  water  as 
bottom  reflectance.  Some  way  must  be  found  to  eliminate  these  po¬ 
tential  erroneous  interpretations  if  Landsat  is  to  be  effective  in 
nautical  charting  applications .  Since  most  of  the  factors  leading  to 
erroneous  interpretations  are  transient,  comparison  of  images  from 
different  dates  can  often  reveal  their  transient  nature  and  permit 
correct  classification  of  the  image  patterns. 

The  Landsat  Sensor  System 

Characteristics  of  the  Landsat  sensor  systems  have  been  well 
described  by  Taranik  (1973) .  The  sensor  of  primary  interest  in  this 
study  was  the  multispectral  scanner  (MSS)  recording  reflected  solar 
energy  in  four  wavelength  bands,  with  the  return- beam  vidicon  sensors 
(RBV)  of  secondary  interest.  The  spectral  bands  associated  with  each 
of  these  sensors  are: 
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Energy  received  at  the  sensor  in  each  of  these  bands  is  tele¬ 
metered  to  ground  receiving  stations.  Data  from  the  four  channels  of 
the  MSS  are  converted  to  digital  values,  each  representing  the 
average  brightness  of  the  scene  within  the  instantaneous- fie Id- of- 
view  (IFOV)  of  the  scanner  in  a  particular  spectral  band,  at  a 
particular  moment  in  time.  These  digital  values  can  be  processed 
in  digital  form  or  converted  to  an  image  form  for  analysis  by  con¬ 
ventional  image  interpretation.  Each  digital  value  represents  an 
area  in  the  scene,  called  a  picture  element  (pixel),  with  an 
assigned  brightness  value  ranging  from  0  through  127  in  Sands  4,  5, 
and  6,  and  from  0  through  63  in  Band  7.  Since  image  products  cannot 
show  more  than  about  20  detectably  different  shades  of  grey,  the  full 
information  content  of  the  record  can  seldom  be  extracted  from  the 
image  form,  alone.  At  the  present  state  of  the  art  in  digital  analysis 
of  Landsat  data  human  interpreters  can  do  a  better  job  of  interpreting 
contextual  information  than  can  computers,  and  joint- human  machine 
interpretations  are  frequently  required  to  maximize  information  ex¬ 
traction.  In  this  study,  human  interpretation  of  multitemporal  data 
was  used  to  identify  image  patterns  associated  with  non-bottom 
features  so  that  these  patterns  would  not  be  erroneously  interpreted 
as  shoals  or  shallow  water  areas.  Computer  processing  of  the  digital 
data  to  obtain  relative  depth  information  was  conducted  in  association 
with  the  human  interpretation  of  the  images.  Computer  processing  was 
performed  at  The  University  of  Michigan  Computer  Center  through  the 
Michigan  Terminal  System  (MTS).  Image  histograms,  statistics,  direct 
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output  of  signal  values  (digital  numbers),  and  line-printer  "grey-maps" 
were  produced  with  the  11  LINE  software  package  developed  at  the 
Environmental  Research  Institute  of  Michigan  (ERIM) .  Statistical 
programs  designed  for  the  Michigan  Interactive  Data  Analysis  System 
(MIDAS)  were  used  to  establish  parameters  for  transforming  geographic 
positions  from  a  Mercator  chart  to  Landsat  column  (line)  row  (point) 
pixel  coordinates,  and  for  other  statistical  operations  required  in 
depth  calculation. 


MULTITEMPORAL  ANALYSES 

When  confronted  with  a  variety  of  images,  the  interpretor  should 
first  look  for  known  features,  or  known  characteristics  before  proceed¬ 
ing  with  the  interpretation  of  unknown  features.  Initial  orientation 
can  be  achieved  by  associating  features  seen  in  Figure  1  with  the  in¬ 
dex  map  of  the  area  (Figure  4).  Figure  4  shows  the  Cape  Cod  peninsula 
and  the  surrounding  bodies  of  water.  The  Atlantic  Ocean  is  to  the 
east,  Cape  Cod  Bay  to  the  north,  Nantucket  Sound  to  the  south  separates 
the  ten- kilometer  wide  peninsula  from  the  two  large  islands  of  Martha's 
Vineyard  and  Nantucket..  Highly  reflective  quartz- sand  dunes,  created 
from  detrital  wind  blown  sand  appear  to  the  north  and  northeast  of 
Provincetown  and  on  Sandy  Neck,  the  spit  which  partly  enclosed  Barn¬ 
stable  Harbor.  Smooth  mature  coastlines  surround  the  Cape  Cod  peninsula 
and  the  associated  islands.  Recurved  spits  such  as  Monorooy  and  Race 
Point  have  grown  from  the  deposition  of  sediment  by  littoral  drift. 
Baymouth  bars  and  spits  have  grown  across  the  mouths  of  estuaries, 
often  sealing  them  off  from  the  open  sea. 
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A  feature  which  appears  in  the  waters  directly  above  Great  Point 
on  Nantucket  Island  is  not  a  spit.  It  is  called  Point  Rip,  is 
subaqueous,  and  is  curved  toward  the  Atlantic  Ocean  not  the  mainland. 
It  marks  the  location  of  the  remnants  of  a  glacial  kame  deposit, 
which  originally  provided  the  sediments  for  Coskata  Beach  and  Great 
Point  (Tisdall  and  El-Baz,  1979). 

A  most  distinctive  feature  is  the  cuspate  shoreline  within  the 
harbor  of  Nantucket.  The  origin  and  mode  of  formation  of  these  cusps 
have  been  disputed  for  some  time.  Farouk  El-Baz  and  Tracey  Tisdall 
(1979)  suggest  the  following  model  as  most  plausible:  (1)  The  shore 
was  modified  by  the  wind  system  particular  to  the  area;  (2)  The 
harbor  is  oriented  parallel  to  the  prevailing  winds  from  the  south¬ 
west;  (3)  The  dominant  storm  wind  direction  is  from  the  northeast; 
and  (4)  It  is  primarily  the  action  of  the  opposing  winds  at  high 
angles  to  the  beach  that  has  produced  the  cuspate  configuration. 
Leonard  and  others  (1976)  felt  that  "longshore  processes  acting  in 
both  directions  erode  sediment  from  the  centers  of  the  concavities 
and  transport  it  to  the  spit  ends  where  it  is  deposited  as  subaqueous 
bars,  the  upwind  spit  preventing  longshore  drift  in  the  center  of 
the  concavity."  This  process  creates  the  symmetry  of  the  cusps  and 
the  lengths  of  the  points  are  determined  by  the  relative  strength  of 
currents.  The  distance  between  the  cusps  is  thought  to  depend  upon 
the  harbor  width.  Equilibrium  of  processes  is  necessary  for  the 
stability  of  the  beach  morphology.  Although  the  beach  has  been 
stable  for  some  time,  there  is  some  evidence  that  an  inlet  formerly 
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existed  at  Haulover  3each,  causing  decreased  tidal  flow  in  the  harbor. 
This  caused  an  imbalance  between  accretion  on  the  spits  and  erosion  by 
the  tidal  currents.  The  spits  then  grew  long  enough  to  reach  the 
opposite  shore.  Now  that  the  inlet  is  closed,  equilibrium  has  again 
been  regained.  The  remnant  bars  crossing  the  harbor  can  be  seen  in 
Figure  5. 

Shoal  Detection 

Light- toned  areas  in  the  waters  around  Cape  Cod  may  represent 
bottom  features  due  to  upwelling  reflectance  through  the  water;  or, 
on  the  other  hand,  they  may  represent  only  temporal  "shoal 
apparitions"  of  phytoplankton  concentrations,  sediment- laden  water , 
clouds  and  haze,  or  other  surface  effects,  such  as  sea  state  and  sun- 
glint.  The  prudent  navigator  would  not  be  in  danger  from  interpre¬ 
tation  errors  of  light- toned  water;  however,  errors  in  interpretation 
of  darker- toned  water  due  to  light  absorbing  vegetation  and  sediments 
could  be  a  real  danger  in  shallow  seas  because  these  factors  can  make 
water  appear  deeper  than  it  really  is.  Several  features  are  con¬ 
sistently  represented  on  multitemporal  images  and  are  also  designated 
as  bathymetric  shallows  on  published  charts. 

Shoal  areas  around  Nantucket  can  be  detected  in  the  waters  to 
the  north  and  to  the  east  of  the  island.  An  enlargement  of  the  17  July 
1974  Landsat  image  of  Cape  Cod  (Figure  6),  overlayed  by  a  portion  of 
the  National  Oceanic  and  Atmospheric  Administration  (NOAA)  National 
Ocean  Survey  (NOS)  bathymetry  chart  No.  13237  (26th  Ed.,  27  January 
1979,  from  1957  survey  data),  illustrates  how  closely  the  lighter- toned 


L8- 


pattems  in  Che  image  correspond  to  Che  bathymetric  highs  of  the 
shoal  areas. 

When  the  attention  of  the  reader  is  focused  upon  Figure  1,  other 
light- toned  areas,  matching  other  areas  of  shallow  bathymetry,  can  be 
detected.  The  gentry  sloping  Billingsgate  Shoal,  for  example,  can  be 
seen  west  of  Jeremy  Point,  a  southern  extension  of  Great  Island  in 
Wellfleet.  The  channel  leading  into  Wellfleet  Harbor  is  also  evident. 
Bars  and  channels  of  ebb  tidal  deltas  can  be  seen  outside  Nantucket 
Harbor,  Barnstable  Harbor,  and  Chatham  Harbor.  Within  Barnstable 
Harbor  are  flood  tidal  deltas.  These  features  are  also  detectable 
but  are  less  apparent  within  Pleasant  Bay.  The  area  of  Brewster  Flats 
east  of  Barnstable  Harbor  and  extending  further  eastward  to  Eastham 
is  also  highly  visible.  The  area  is  rather  regular  in  appearance  when 
contrasted  to  the  sinuous  patterns  of  shoals  around  Nantucket  and 
within  Nantucket  Sound  north  of  Tuckernuck  and  Muskeget  Islands. 

This  is  a  result  of  the  lower  energy  of  wave  turbulence  due  to 
shorter  wind  fetch  within  Cape  Cod  Bay  compared  to  the  coast  towards 
the  Atlantic  Ocean.  To  the  south  of,  and  adjoining  Nantucket  Island 
at  its  midsection,  is  another  remnant  of  a  glacial  kame  called  Miacomet 
Rip.  Within  Nantucket  Sound,  on  the  west  side  of  Monomoy  Island,  is 
the  Commons,  an  oLd  tidal  delta.  This  bathymetric  feature  will  be 
discussed  from  another  point  of  view  later  in  this  thesis. 

The  offshore  bar  called  Peaked  Hill  3ar  can  be  detected  north  of 
Provincetown  with  aid  of  magnification.  The  resolution  of  the  Landsat 
MSS  (79  meters)  is  simply  not  adequate  to  discern  smaller  offshore 
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bars  without  additional  knowledge  and  ancillary  data.  It  can  be 
shown  that  smaller  features  in  the  water  can  be  detected  with  digital 
data,  but  the  features  may  be  just  too  difficult  to  recognize.  The 
narrow-band  multispectral  quality  of  Landsat  imagery  was  found  to  be 
of  significant  benefit  over  broad-band  panchromatic  imaging  systems. 

The  higher  resolution  of  Landsat- 3  RBV  (40  meters)  images  have 
significant  advantage  over  the  MSS  for  recognition  of  submarine 
features.  By  comparing  Figure  1  to  Figure  5  it  is  evident  that  much 
greater  detail  in  submarine  morphology  can  be  more  easily  delineated. 
In  the  Landsat- 3  RBV  image,  breaking  waves  can  be  seen  approximately 
over  the  location  of  the  shoals.  Offshore  bars  are  exposed  along 
the  southern  shore  of  Nantucket.  Much  more  detail  of  under  water 
channels  and  shoals  can  be  seen  near  the  islands  of  Tuckernuck  and 
Muskeget,  west  of  Nantucket.  Features  appear  in  the  Commons  which 
may  be  parallel  bars  oblique  to  the  shoreline  of  Monomoy.  The 
tidal  flats  parallel  to  the  shoreline  north  of  Nantucket  appear 
cross-hatched  and  may  indicate  the  presence  of  bars,  although  these 
features  do  not  appear  on  the  latest  edition  of  NGAA's  (NOS)  chart 
No.  13237  (26th  Ed.,  27  January  1979). 

The  resolution  of  the  S-190B  Earth  Terrain  Camera  (ETC)  "high- 
resolution"  color  photograph  (SL- 3-86- 314)  in  Figure  3  is  21  meters 
at  a  scale  of  1:250,000  (enlarged  from  an  original  scale  of  1:950,000). 
The  photographs  taken  over  Cape  Cod  are  superior  in  clarity  and  reso¬ 
lution  to  Landsat  MSS  or  Landsat- 3  RBV  imagery.  Ship  wakes,  exposed 

bars,  breaking  waves,  and  even  subaqueous  sand  waves  lying  above  the 
Nantucket  Shoals  are  clearly  visible  in  the  original  transparency. 


Sea  State 


The  extreme  variety  of  patterns  of  light  and  dark  tones  are 
represented  by  multi-date  band  4  Landsat  MSS  images  in  Figures  7-10. 
Careful  analysis  of  the  patterns  to  the  east  of  Nantucket  revealed 
a  clue  that  the  patterns  were  related  to  sea  state.  The  edges  of  the 
pattern  in  the  28  May  1975  Landsat  image  (Figure  7)  are  very  different 
from  the  17  July  1974  Landsat  image  (Figure  1).  In  Figure  1  the 
eastern  edge  is  sharp  and  the  western  edge  is  diffuse.  The  edges  of 
the  same  pattern  in  Figure  7  are  just  the  opposite,  indicating  a 
change  of  wind  direction  upon  the  sea  surface.  A  dramatic  change  is 
apparent  in  the  9  March  1978  image  (Figure  9),  where  distinct  patterns 
of  tone  are  barely  detectable.  The  wind  direction  was  from  the 
north- northeast  (20°)  for  the  17  July  1974  image;  north- northwest 
(330°)  for  the  28  May  1975  image;  and  no  wind  was  present  at  Nantucket 
Tower  Weather  Station  at  the  time  of  the  9  March  1978  image.  Other 
patterns  in  Figures  7-10  generally  seem  to  correlate  to  charted 
bathymetric  highs,  with  variations  due  to  wind  direction,  wind  speed, 
and  wave  swell.  Sea  state  conditions  could  cause  sediment  to  be 
entrained  directly  above  the  higher  points  of  the  shoals.  The  amount 
of  sediment  suspended  is  a  function  of  wave  height  to  water  depth. 
Therefore,  the  height  of  the  tide  is  also  important. 

The  9  March  1978  MSS  image  is  exceptional  for  many  reasons.  It 
is  one  of  only  a  few,  cloud- free,  wind- free  images  obtained  by  Landsat 
over  Gape  Cod  in  the  five-year  sampling  period.  The  absence  of  wind 
generated  surface  effects  leaves  the  water  areas  in  the  image  devoid 
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of  the  striking  multitude  of  patterns  that  characterize  the  other 
images;  however,  some  patterns  are  visible.  They  are  the  shoals, 
whose  upwelling  reflectance  is  not  masked  by  water  turbulence.  They 
are  less  prominent  because  proportionally  less  direct  sunlight 
illuminates  and  reflects  from  them  to  the  sensor.  This  is  due  to 
the  low  sun  elevation  (above  the  horizon)  and  increased  path  radiance, 
than  would  be  the  case  for  images  with  higher  sun  elevations.  The 
light  reaching  the  bottom  features  and  subsequently  reflected  is  also 
attenuated  by  the  level  of  water  (2.4  meters)  above  mean  low  tide. 

The  reader  will  note  that  the  mainland  and  Martha's  Vineyard  are 
much  brighter  than  Nantucket.  The  difference  in  brightness  is  not  a 
result  of  special  processing,  or  of  radiometric  discrepancies,  because 
the  Landsat-3  RBV  images  from  the  same  day  and  time  also  contain  the 
striking  difference  in  brightness  between  the  land  surfaces.  3oth  the 
MSS  and  RBV  images  were  taken  in  late  winter,  when  Nantucket  had 
apparently  escaped  a  snow  fall  which  covered  the  mainland  and  Martha's 
Vineyard. 

Radiometric  discrepancies  do  occur  in  the  Landsat-2  images 
created  in  the  period  between  22  January  to  16  July  1975  (Figures  7 
and  10).  The  preflight  calibration  data  used  did  not  fully  utilize 
the  dynamic  range  capabilities  of  the  sensor  and,  therefore,  had  a 
different  response  compared  to  other  Landsat  responses  (Colvocoresses , 
1977). 

The  images  appear  lighter  and  do  not  have  as  much  contrast.  They 
were  electronically  "clipped;"  e.g.,  the  digital  number  values  were 


adjusted  by  about  three  digital  numbers  at  the  high  end  and  one  digital 
number  at  the  low  end  of  the  128  discrete  grey- levels  for  MSS  3ands  4, 
5,  and  6  and  64  discrete  grey- levels  for  MSS  Band  7.  This  was  done 
partly  to  reduce  striping  caused  by  differences  in  the  output  of  the 
individual  detectors. 

Blizzard  and  Landfora  Changes 

Another  hypothesis  that  was  considered  in  the  interpretation  of 
the  9  March  1980  image  was  the  possibility  that  the  shoal  formations 
might  have  been  altered  and  made  less  visible  by  the  catastrophic 
blizzards  of  1978.  Damage  from  that  series  of  blizzards  (especially 
the  6-7  February  1978  storm)  has  been  estimated  at  over  one  billion 
dollars.  Piers,  jetties,  breakwaters,  and  beach- front  homes  were 
destroyed  by  the  nine  meter  waves,  generated  by  50  knot  winds  that 
blew  for  up  to  33  hours. 

The  storm  severely  impacted  the  coastal  areas  from  Rhode  Island 
to  Maine.  On  York  Beach,  Maine,  an  ancient  wreck  was  uncovered  by 
the  storm- -reminiscent  of  the  story  of  the  Sparrow-Hawk.  The  exposed 
hull  of  the  small  ship  was  fifty  feet  long  and  between  300-600  (sic.) 
years  old  (Phaneuf,  1978).  At  Provincetown,  the  sea  separated  the 
northern  tip  of  Cape  Cod  from  Herring  Cove  Beach.  Some  residents  lost 
as  much  as  600  feet  of  shoreline  to  the  ocean  (Phaneuf,  1978).  The 
fishing  fleet  of  Scituate  was  stranded  and  a  tanker  ship  was  grounded 
at  Marblehead  (Anonymous,  1978).  Monomoy  Island,  which  had  been  a 
single  island  for  many  years  was  breached  by  the  sea  and  is  now  sepa¬ 
rated  into  two  islands  (Figures  11  and  12).  The  Commons,  the  tidal 
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delta  mentioned  earlier,  is  a  result  of  a  break  in  Monomoy  Island 
approximately  150  years  ago  (Tisdall  and  El-Baz,  1979).  The 
separation  of  Monomoy  Island  from  the  mainland  is  reported  to  have 
occurred  in  1960.  The  separation  was  triggered  by  the  construction 
of  a  causeway  between  Morris  Island  and  Chatham  proper  (Tisdall  and 
El-Baz,  1979).  A  new  tidal  delta  has  already  begun  to  form.  It  is 
more  complex  than  the  Commons  delta,  because  it  has  both  an  ebb  and 
a  flood  component. 

Littoral  Drift 

Significant  changes  were  noticed  in  the  subaerial  land  forms 
that  were  not  recorded  on  the  charts.  Several  spits  have  an  altered 
shape  or  join  small  islands.  Jeremy  Point  has  now  joined  Billingsgate 
Island  by  extending  its  length  500  meters  since  1967  (Tisdall  and 
El-Baz,  1979).  These  changes  were  not  due  to  change  in  tide  level 
on  different  Landsat  images.  Nauset  Spit  migrates  southward  at  a 
rate  of  about  90  meters  per  year  (Strahler,  1966).  Sediment  has  been 
removed  from  the  midsection  of  Monomoy  by  wave  refraction  and  trans¬ 
ported  southwards  by  littoral  drift  resulting  in  an  increase  in  its 
length  (Williams,  1979;  Odale,  et  al.,  1971).  The  growth  can  be  seen 
in  the  trend  of  sequential  arcs  which  mark  past  changes  in  the  de¬ 
position  process.  Another  indicator  of  current  direction  is  the 
accumulation  of  sediment  on  the  updrift  side  of  jetties,  groins,  and 
similar  barriers  as  seen  on  (Figure  12)  the  northern  boundary  of 
Nantucket  Sound  and  verified  by  Strahler' s  (1966)  published  longshore 
drift  current  map  (Figure  13) .  Sediment  transport  in  the  area  of  Cape 
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Cod  is  primarily  accomplished  by  littoral  drift,  except  in  the  dune 
areas  at  Provincetown  and  at  Sandy  Neck. 

Turbidity 

Turbidity  is  defined  as  the  suspension  in  the  water  of  fine  silt 
sized  minerals  biological  particulate  matter,  and  "yellow  substance" 
(from  the  German,  "gelbstoffe") ,  (a  strong  blue- absorbing  complex  of 
decomposing  vegetative  matter  that  resembles  tannic  acid).  Better 
bathymetric  data  can  be  obtained  if  areas  of  suspected  turbidity  can 
be  detected  and  then  edited  for  alternative  processing  by  use  of 
more  appropriate  attenuation  coefficients.  It  is  difficult  to  visually 
detect  changes  of  water  turbidity,  but  general  guidelines  can  help 
Landsat  users  evaluate  areas  where  water  clarity  change  is  suspected. 

Very  fine-grain  sediment  can  remain  suspended  for  extremely  long 
periods  of  time.  "Yellow  substance"  is  a  complex  product  of  decay, 
and  is  more  prevalent  in  waters  with  a  land  source.  In  low  wave 
energy  environments,  such  as  the  Amazon  continental  margin,  fine- 
grain  sediment  deposition  can  occur  near  the  source.  A  similar 
deposition  would  generally  occur  in  an  estuarine  environment.  In 
high  wave  energy  areas,  such  as  around  Cape  Cod,  fine- grain  sediment 
is  transported  further  from  the  coast  to  the  continental  shelf  or  deep 
ocean  before  deposition  occurs. 

The  highly  turbulent  waters  caused  by  strong  currents  and  breakers 
over  the  shoals  usually  carries  entrained  sediment.  However,  the  sedi¬ 
ments  of  this  relict  beach  have  been  reworked  for  two  millenia  (6000 
years),  and  are  relatively  coarse  (Rea,  1979).  Sediments  measured  by 
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Stewart  and  Jordan  (1964)  on  Georges  3ank  measured  0.68  milimeters 
in  mean-grain  size  on  Che  ridges.  Recenc  studies  of  suspended  matter 
in  the  vicinity  of  Nantucket  inspired  by  the  Argo  Merchant  oil  spill, 
indicate  that  in  surface  waters  the  particulate  count  is  1  miligram/ 
liter  and  at  the  bottom  2  miligram/liter .  After  severe  winter  storms 
under  conditions  of  high  sea  state  in  the  area  south  of  Nantucket 
Island,  the  surface  waters  contained  3  miligram/liter,  the  bottom 
15  miligram/liter.  Non-mir.eral  material  is  combustible,  therefore, 
to  determine  how  much  particulate  matter  is  biological  in  origin, 
samples  are  burned  and  weighed.  Combustibles  represent  50  percent 
of  the  total  particulate  weight.  Turbidity  was  also  measured  with  a 
cransmissometer  and  was  found  to  be  negligible  over  the  shoals  (Bothner, 
et  al.,  1980). 

It  is  reasonable  to  expect  that  heavy  particles  would  settle  out 
as  a  minimum  threshold  of  turbulence  is  reached.  Sediment  transport 
occurs  within  one  meter  of  the  shoal  surface  (Stewart  and  Jordan,  1964). 
Local  fishermen  and  boat  operators  report  that  the  water  is  clear 
enough  to  see  the  bottom  (Ryder,  1979). 

Sediment  plumes  can  usually  be  found  in  the  Landsat  images  near 
the  mouths  of  estuaries,  rivers,  channels,  and  along  spits  where  areas 
of  turbulence  are  occurring.  A  small  plume  is  visible  at  the  northern 
end  of  Cape  Cod  Canal  (Figure  14) .  Sediment  plumes  occur  at  the  point 
where  longshore  currents  moving  in  opposite  directions  converge,  such 
as  midway  along  the  southern  coast  of  Martha's  Vineyard  (Figure  15)  and 
midway  along  the  Atlantic  coast  of  Cape  Cod  (Figure  15) . 
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The  contrail  seen  in  Figure  16  is  a  relative  measure  of  the 
turbidity  in  the  water.  Contrail  shadows  would  not  be  visible  over 
deep,  clear  water  (theoretically),  because  only  skylight  is  reflected, 
and  direct  solar  irradiation  is  completely  attenuated  in  deep  water. 

The  image  suggests,  therefore,  that  there  is  either  turbidity  present 
in  the  water  to  reflect  direct  solar  illumination  from  the  suspended 
particles  or  that  there  is  a  surface  effect.  The  shadow  represents 
the  difference  between  strong  solar  illumination  which  has  been 
interrupted  and  skylight  illumination  which  is  always  present. 

Biological  Factors 

The  fishing  grounds  off  Cape  Cod  and  Georges  Bank  are  some  of  the 
most  productive  in  the  world.  Plankton  is  the  primary  link  in  the  food 
chain.  Plankton  are  not  distributed  evenly,  rather  they  appear  on  the 
surface  in  patches  with  concentrations  accounting  for  about  10  percent 
of  the  particulate  matter.  Plankton  blooms  occur  most  commonly  in 
March  and  September  oftentimes  causing  substantial  fish  kills  (3ush, 
1979).  No  plankton  communities  could  be  detected  with  certainty  in 
the  images,  although  the  water  mass  boundary  seen  in  Figure  10  may  be 
of  biological  origin.  A  further  discussion  of  the  detection  of  phyto¬ 
plankton's  chlorophyll- A  will  follow  later  in  the  paper. 

No  marine  vegetation  could  be  detected  in  the  water,  although  it 
may  be  present.  Typically,  there  is  none  visible  on  or  near  the  sur¬ 
face  (such  as  kelp),  although  marine  plants  foul  bottom  planted 
scientific  instruments  with  up  to  six  inches  of  growth  in  less  than 
three  years  (Bothner,  et  al.,  1980).  Vegetation  in  the  channels  between 
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che  shallower  shoal  areas  is  unlikely  because  of  the  strong  currents 
of  the  rips. 


Sea  Foam 

Sea  foam  is  a  factor  that  should  be  considered  in  the  interpre¬ 
tation  of  the  imagery.  Quantitative  data  of  the  coverage  and  thickness 
of  foam  does  not  exist.  Its  distribution  is  dependent  upon  the  water 
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chemistry,  wind  and  wave  interaction,  spray,  and  bathymetry.  Spilling 
waves  seem  to  generate  more  foam  than  breaking  waves.  Sea  foam  is 
observed  at  the  breakers  and  is  common  to  shallows  (Limeburner,  1979), 
where  it  may  be  extensive  depending  upon  the  wind  speed.  The  per¬ 
centage  of  area  covered  has  been  roughly  estimated  to  be:  0  percent 
up  to  15  knot  winds;  10  percent  at  15  knot  or  higher;  40  percent  at 
20  knot  winds  by  Richard  Limeburner  (1979).  This  estimate  may  be  very 
high  at  the  top  end  but  it  is  significant  that  up  to  15  knot  winds 
that  there  is  little  foam  observable.  This  appears  to  agree  with 
photographs  and  sea  state  description  data  examples  in  the  3owditch 
(U.S.  Defense  Mapping  Agency  Hydrographic/Topographic  Center).  It 
also  shows  agreement  with  low  altitude  aerial  photographs  of  the 
stranded  Argo  Merchant  on  Fishers  Rip  (Milgram,  1977)  where  few  white- 
caps  and  foam  are  seen  even  at  wind  speeds  greater  than  20  knots . 

An  interesting  phenomenon  is  that  the  observed  major  pattern 

which  is  seen  over  the  shoals  appears  in  all  four  spectral  bands  of 

Landsat.  Figure  18  of  17  July  1974  shows  the  patterns  in  MSS  Band  7 

(800-1100  nanometers).  Landsat  MSS  Band  7  has  such  a  high  attenuation 
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coefficient  that  its  water  penetration  is  negligible  (4.66  meters) 
(Wolf  and  Zissis,  1978).  In  addition,  Seasat  Synthetic  Aperture 
Radar  (SAR)  backscatter  images  (Figure  19)  show  patterns  over  the 
shoals  that  are  very  similar  to  those  from  Landsat  (Figure  18) .  The 
L-band  wavelength  (23  centimeters)  of  the  Seasat  Radar  does  not  pene¬ 
trate  the  water.  Therefore,  patterns  observed  must  be  due  to  surface 
effects,  perhaps  a  modulation  of  the  surface  wave  structure,  rather 


than  (or  in  addition  to)  turbidity  within  the  water.  This  phenomenon 
appearing  over  the  shoals  may  be  very  promising  for  charting  if 
shallow  bathymetry  is  directly  linked  to  observable  surface  effects. 
Hazards  to  shipping  could  be  detected  in  spite  of  turbidity.  It  may 
be  of  great  value  for  navigation  of  muddy  rivers.  These  radar 
applications  are  currently  being  investigated  by  Robert  Shuchman  at 
ERIM  (1980,  pers.  comm.). 

Surface  Manifestations  of  Internal  Waves 

Perhaps  one  of  the  more  remarkable,  more  easily  recognized,  yet 
not  absolutely  understood  phenomenon  is  seen  in  Figure  IB.  The  long 
periodic  wavelike  features  have  been  observed  on  spacecraft  imagery, 
including  Landsats  1  and  2,  Defense  Meterological  Satellite  Program 
(DMSP)  satellite;  National  Oceanic  and  Atmospheric  Administration 
(NOAA)  weather  satellites,  Skylab;  Apollo-Soyuz;  and  high  flying  U- 2 
aircraft,  since  1972.  Such  features  usually  appear  on  continental 
shelves  of  the  world  and  are  generally  oriented  parallel  to  the  con¬ 
tinental  slope,  although  sometimes  with  variations  which  are  apparently 
influenced  by  topographic  and  bathymetric  features  (Sawyer  and  Apel, 
1976).  The  circular  feature  on  Figure  1  is  one  result  of  such  in¬ 
fluence.  The  morphology  below  this  visible  phenomenon  is  a  basin 
which  has  one  prominence  on  its  southern  edge  as  indicated  on  Figure 
30  of  the  following  section.  They  have  also  been  observed  in  deep 
ocean  areas  (Apel,  et  al.,  1976)  and  off  the  North  .American,  African 
East  and  West  coasts,  in  the  Baltic  Sea,  in  the  Suiu  Seas,  in  the  Gulf 
of  Mexico,  and  in  the  Caribbean  Sea  (Apel,  et  al.,  1975). 


They  are  theorized  to  be  the  visible  surface  manifestations  of 
internal  gravity  waves.  These  internal  wave  structures  result  from 
stratification  within  the  water  body  of  two  different  densities  of 


water.  The  density  difference  can  result  from  variations  in  tempera¬ 
ture  and/or  salinity.  These  surface  manifestations  are  distinguished 
by  cheir  large  scale  (Apel,  et  al.,  1976)  and  are  most  frequently 
seen  during  the  summer  months . 

Observed  from  the  surface  they  appear  as  long,  narrow 
quasi-parallel,  smooth  streaks  on  the  water  when  the 
wind  is  light.  "They  move  slowly  toward  the  coastline, 
one  after  the  other,  each  locating  below  it  the  point 
on  the  internal  wave  halfway  between  a  crest  and  the 
following  trough,  as  this  is  the  region  where  the  cur¬ 
rents  at  the  surface,  due  to  the  internal  wave  below, 
converge.  The  natural  oils  damp  the  capillary  waves 
and  produce  a  smooth,  slick  appearance."  (Neumann  and 
Pierson,  1966). 

This  action  is  thought  to  increase  specular  reflection  and  to  de¬ 
crease  the  diffuse  scattering  over  the  convergence  region.  "A  second 
theory  predicts  that  the  small  waves  are  concentrated  in  the  con¬ 
vergence  regions  due  to  the  wave-current  stresses,  thereby  decreasing 
the  specular  reflection  and  increasing  the  diffuse  scattering  over 
the  convergence  region"  (Apel,  et  al.,  1976). 

The  qualitative  characteristics  of  the  surface  manifestations  of 
internal  waves  when  studied  on  spacecraft  images  are  summarized  by 
Apel,  et  al. ,  (1976) ; 

(1)  The  waves  occur  in  groups  or  packets  of  3  to  5  kilometers 
usually  landward  of  the  continental  slope  and  separated  by  dis¬ 
tances  which  are  on  the  order  of  either  15  or  30  kilometers; 
taken  together  with  both  the  observed  and  calculated  phase 
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velocities,  these  facts  suggest  a  semidiurnal  or  a  diurnal 
origin; 

(2)  The  crests  are  nearly  always  oriented  parallel  to  the 
bottom  topography  or  can  be  loosely  associated  with  some  topo¬ 
graphic  feature  seaward  of  their  observed  position,  or  both; 

(3)  The  wavelengths  fall  between  200  and  10,000  meters,  depend¬ 
ing  on  the  geographical  area;  within  a  given  packet,  there  is  a 
monotone  decrease  in  wavelength,  varying  from  X  at  the  front  to 

X  at  the  back  of  the  group ; 

o 

(4)  The  lengths  of  the  crests  fall  between  a  very  few  and  per¬ 
haps  100  kilometers  with  a  decrease  in  crest  length  occurring 
from  front  to  back  of  the  group; 

(5)  The  widths  of  the  slicks  are  often  small  compared  to  the  lead 
wavelength; 

(6)  The  crests  are  curved  in  a  horizontal  plane  with  their  con¬ 
vex  sides  pointed  in  the  direction  of  propogation;  the  radii  of 
curvature  range  from  essentially  infinity  to  a  few  kilometers. 

(7)  As  the  packet  progresses  up  on  the  continental  shelf,  there 
is  some  evidence  of  a  continued  increase  in  the  wavelengths 
throughout;  an  accounting  for  this  may  be  had  by  a  combination 
of  linear  dispersion  and  nonlinear  effects  akin  to  solitary  wave 
behavior. 

Unfortunately,  simultaneous  identification  of  internal  waves  on 
spacecraft  images  and  from  on-site  ship  observation  has  not  yet  been 
completely  successful;  however,  all  data  gathered  from  shipboard 
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thermistors  and  sonar  in  the  New  York- to- 3ermuda  Remote  Sensing 

Experiment  (NYBERSEX)  by  Apel,  et  al.,  (1976)  are  consistent  with 

the  characteristics  as  seen  in  associated  imagery. 

Internal  waves  can  affect  ships.  The  density  discontinuity 
between  water  strata  can  be  quite  sharp  under  certain  cir¬ 
cumstances.  "An  internal  wave  motion  at  the  interface  is 
easily  excited  as  very  little  energy  is  needed  to  cause  a 
large  internal  oscillation  that  has  at  the  same  time,  very 
little  kinetic  and  potential  energy.  The  currents  are  slow 
and  work  done  in  displacing  the  denser  mass  of  water  is 
small  as  it  displaces  only  a  slightly  less  dense  mass  of 
water.  If  the  density  discontinuity  is  near  the  keel  of 
a  ship,  it  is  observed  chat  large  amounts  of  power  are  needed 
to  get  the  ship  underway  and  to  maintain  speed  if  the  ship 
is  going  slowly.  This  is  the  phenomenon  of  'dead  water,' 

For  the  same  power,  the  ship  moves  more  slowly  as,  in 
addition  to  making  bow  waves  and  a  wake  and  overcoming 
friction,  she  must  also  generate  internal  waves"  (Neumann 
and  Pierson,  1966  .  • 

Internal  waves  pose  no  threat  to  the  safety  of  surface  vessels, 
but  are  considered  a  potential  hazard  to  submarines.  Submarines 
maintain  their  underwater  vertical  position  by  ballasting.  If  a 
submerged  vessel  over- ballasted  upon  encountering  an  upward  accelera¬ 
tion  due  to  an  internal  wave,  and  was  subsequently  caught  in  the 
downward  acceleration  phase  of  the  wave  cycle,  it  could  conceivably 
exceed  the  ship's  maximum  depth  limit  or  strike  bottom  before  com¬ 
pensation  of  the  over- ballast  could  be  accomplished.  This  scenario 
has  been  suggested  as  a  possible  cause  of  the  sinking  of  the  American 
submarine  Scorpion  and  the  subsequent  loss  of  all  her  crew  (Naval 
Environmental  Prediction  Research  Facility,  1979  >. 


Tidal  Overfalls  on  the  Shoals 


Another  factor  that  should  be  tested  for  detectability  by  multi¬ 
temporal  Landsat  data  is  the  phenomenon  of  water  motion  variously  called 
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"intermictent  tide  rip,"  "a  moving  wave,"  and  "a  confluence  of  ocean 
currents,"  which  has  been  reported  as  a  pecularity  of  the  area  by 
Stewart  and  Jordan  (1964).  Stewart  and  Jordan  (1964)  found  that 
there  is  a  visible  tidal  overfall  that  occurs  when  a  large  volume 
of  fast  moving  water  when  it  is  forced  upward  and  over  the  crests  of 
the  sand  shoals  when  the  tidal  current  flow  direction  is  perpendicular 
to  the  orientation  of  the  shoal,  as  at  the  time  of  the  Landsat  image 
of  17  July  1974.  Perhaps  this  phenomenon  of  water  motion,  which  has 
such  a  striking  visual  impact  on  ground  observers  (compared  to  the 
impact  on  surface  observers  of  manifestations  of  internal  waves), 
might  contribute  to  the  patterns  seen  above  the  shoals  on  Landsat  and 
other  satellite  images. 

Pollock  Rip  Channel  was  chosen  as  the  site  of  this  investigation, 
because  tidal  current  information  was  readily  available  and  highly 
accurate,  because  it  is  the  location  of  a  NOAA  observation  station, 
whereas  Nantucket  Shoals  has  no  stations  which  could  provide  direct 
observations . 

Table  4  lists  the  results  of  the  study  of  image  prints,  currents, 
and  winds.  This  part  of  the  study  was  inconclusive  in  regard  to  the 
detection  of  tidal  overfall,  because  the  factors  of  windspeed,  wind 
direction,  and  radiometric  rendition  of  the  printed  photographs  far 
outweighed  any  detectable  visual  contribution  to  the  brightness  of 
the  patterns.  Fifteen  knot  winds  from  the  north  and  northeast  gener¬ 
ated  the  brightest  patterns  over  Pollock  Rip  than  any  of  the  other 
images,  with  the  exception  of  one  influenced  by  errors  of  image  proces¬ 


sing  . 
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Future  studies  should  be  pursued,  especially  with  regard  to  SAR 
radar  imagery,  because  of  the  documented  alteration  of  the  surface 
wave  structure.  A  larger  sample  of  digital  data  should  also  be 
analyzed  and  correlated  with  tidal  current  flow  direction  and  with 
wind  velocity  and  direction,  so  that  a  quantitative  relationship 
might  be  determined. 

Sung lint 

Interpretation  of  what  might  be  the  cause  of  the  lighter  grey- 
toned  water  areas  (Figure  1C)  in  the  upper  right  portion  of  the  image 
as  contrasted  with  the  dark  areas  (Figure  1)  proved  to  be  difficult 
with  the  available  multitemporal  data.  Changes  detected  could  not 
conclusively  be  linked  positively  or  negatively  to  the  three  most 
likely  variables  of:  (1)  white  caps  due  to  sea  state,  (2)  haze,  or 
(3)  sunglint,  when  comparing  the  same  nominal  Landsat  scene. 

Haze  as  a  contributor  to  the  "grey-water"  areas  could  be  safely 
eliminated  on  the  17  July  1974  image  under  study,  because  details  of 
the  internal  wave  manifestations  could  be  seen  clearly.  Details  of 
features  on  land  also  could  be  seen  clearly;  and  the  absence  of  fog 
or  haze  for  land  areas  is  supported  by  weather  logs  which  show  clear 
visibility. 

Sea  state  due  to  wind  waves  and  swell  was  shown  to  be  higher  at 
Chatham  Coast  Guard  Station,  where  areas  of  "grey-water"  are  visible 
in  the  scene;  then  at  the  "dark  water"  location  of  Nantucket  Light¬ 
ship,  where  calm  was  reported  (Table  2).  It  is  probable  that  the  grey 
areas  in  the  upper  right  portion  of  the  image  also  have  rougher  waters. 
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Sunglint  is  a  result  of  specular  reflection  of  direct  solar 
irradiance  upon  the  water.  It  occurs  when  the  angle  of  incidence  is 
equal  to  the  angle  of  view  of  the  sensor.  Sunglint  in  photometric 
photographs  appears  as  a  round  blob  of  light  with  a  bright  core  (Cox 
and  Munk,  1954).  The  pattern  spreads  over  more  area  when  the  sea 
state  is  higher  (Naval  Environmental  Prediction  Research  Facility, 
1979).  Reflectance  is  a  function  of  the  index  of  refraction,  the 
viewing  angle,  and,  because  wind  affects  the  statistics  of  wave 
slopes  on  the  ocean  surface,  it  is  also  a  function  of  windspeed 
(Duntley,  et  al.,  1974).  Sunglint  pattern  on  optical-mechanical 
scanner  images  is  usually  seen  as  a  lighter- toned  strip  the  length 
of  the  image  on  the  side  in  the  direction  of  the  sun;  good  examples 
of  which  can  be  found  in  the  manuals  of  the  Defense  Meterological 
Satellite  Program  (Naval  Environmental  Prediction  Research  Facility, 
1979). 

On  Landsat  images,  because  of  the  satellite's  selective  view 
angles  and  sun  synchronous  orbit,  sunglint  is  minimized.  Therefore, 
sunglint' s  likelihood  as  a  causative  agent  of  the  "grey-water"  patterns 
was  not  obvious,  until  sequential  frames  of  two  photographs  (panchro¬ 
matic;  0.5-0. 6  micrometers)  from  the  S190-A  camera  of  Skylab-3,  taken 
on  12  September  1973,  were  reviewed  (Figures  20  and  21).  The  six- 
lens  multispectral  camera  produced  photographs  with  an  original  scale 
of  1:2,850,000.  These  figures  are  enlarged  to  a  scale  of  1:500,000. 

The  sequence  illustrates  how  the  pattern  changes  with  the  angle  of 
view.  This  change  would  not  occur  with  a  Lambertian  reflector,  such 
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as  haze  or  sea  foam.  The  effects  of  sunglint  should  therefore,  be 
conclusively  determined  by  the  change  of  view  angle  as  the  sensor 
platform  traverses  its  orbital  path;  however,  with  sequential  Landsat- 1 
scenes  (Figure  22-24)  the  anticipated  evidence  of  sunglint  pattern 
did  not  appear.  3ecause  Landsat  is  an  optical-mechanical  scanner 
system,  it  does  not  take  an  instantaneous  image  as  does  a  camera 
system.  The  imaged  pixel  areas  are  in  a  state  of  flux,  and  each 
water  surface  facet  variously  reflects  specularly  or  "twinkles  and 
glitters."  The  rise  in  signal  received,  therefore,  may  be  attributed 
to  either  specular  reflection  or  to  sea  state.  The  range  of  sun 
azimuths  of  the  images  was  not  sufficient  to  notice  any  difference 
in  pattern.  The  RBV  sensor  which  does  capture  an  instantaneous  image 
was  not  used  for  this  same  sequence  of  scene  coverage. 

Summary  of  Multitemporal  Study 

This  portion  of  the  paper  illustrates  the  value  of  interpretations 
of  multi temporal  Landsat  data.  Highly  detailed  information  can  be 
gathered  from  orbital  space  sensors  about  offshore  features,  some  of 
which  do  not  appear  on  charts.  These  results  indicate  that  rapid 
change  in  coastal  and  submarine  landforms  can  be  monitored  by  satellite. 
Catastrophic  changes,  such  as  those  caused  by  a  "100- year  storm,"  can 
be  assessed  quickly,  and  areas  of  survey  priority  can  be  determined. 
Trends  of  regional  changes  can  also  be  predicted.  Space  imagery  is 
particularly  well  suited  to  the  monitoring  of  coastal  environments 
because  of  the  large  areal  coverage  which  can  be  captured  by  a  single 
frame.  Space  imagery  can  amplify  studies  done  during  periods  of 
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limited  areal  coverage  from  either  photographs  taken  by  low  flying 
survey  aircraft  or  from  observations  from  surface  ships. 

A  great  deal  of  information  has  been  obtained  from  evaluation 
of  surface  manifestations  of  internal  waves  on  Landsat  imagery,  some 
of  which  was  not  previously  known.  Significant  information  has  been 
acquired  from  multisensor  data.  The  patterns  seen  over  the  shoals 
are  created  by  a  significant  contribution  from  surface  effects,  such 
as  sea  state,  which  are  detectable  by  Seasat  radar  as  well  as  by  the 
sensors  on  Landsat.  The  hypothesis  that  the  patterns  seen  in  the 
Landsat  image  of  17  July  1974  are  caused  by  "muddy  water"  is  unlikely 
to  be  correct,  because:  (1)  The  shoals  are  composed  of  coarse  sedi¬ 
ments  which  were  of  glacial  outwash  plain  deposition;  (2)  The  sedi¬ 
ments  were  sorted  by  glacial  runoff  water;  (3)  The  sediments  were 
reworked  by  several  subsequent  periods  of  wave  action  when  the  sea 
level  rose  and  fell;  and  (4)  Because  the  fine  grained  sediments  were 
transported  by  wave  and  current  action  to  deeper,  calmer,  water  for 
deposition.  Sunglint  is  not  easily  detected  from  analysis  of  single 
nominal  Landsat  scenes  but  detection  of  sunglint  improves  when 
sequential  scenes  are  studied.  The  amount  of  information  gathered 
without  aid  of  surface  observations  is  significant. 


MULTI SPECTRAL  ANALYSES 


The  human  interpretations  of  remote  sensor  images  described  in  the 
preceding  section  provided  important  insights  into  hydrographic  features 
of  the  Cape  Cod  area.  Edge  definition  of  some  image  patterns  indicates 
breaking  waves  over  submarine  features.  Some  image  patterns  were  found 
to  correlate,  at  least  approximately,  with  charted  bathymetry,  while 
others  were  apparently  due  to  sea  state.  Sea  state  was  found  to  be  a 
contributing  factor  to  the  phenomenon  of  "grey-water”  areas  (Figure  1C), 
and  also  to  the  visibility  of  the  surface  manifestation  of  internal 
waves  (Figure  IB) . 

Based  on  these  human  interpretations  it  appears  clear  that  the 
Landsat  MSS  sensor  can  be  used  to  assist  in  locating  possible  shoal 
areas,  and  to  determine  if  the  observed  phenomena  are  due  to  bottom  re¬ 
flectance,  turbidity,  or  surface  effects.  For  hydrographic  charting 
purposes,  it  is  important  to  know  when  patterns  observed  in  the  image 
are  due  to  bottom  reflectance  or  other  factors.  When  an  observed  pat¬ 
tern  is  not  related  to  water  depth,  it  is  important  to  determine  what 
other  factor (s)  could  account  for  the  image  appearance. 

Since  sea  water  has  negligible  water  penetration  in  the  infrared 
portion  of  the  spectrum  (Wolfe  and  Zissis,  1978),  patterns  seen  in  all 
four  Landsat  bands  should  not  be  due  to  bottom  reflectance.  Image 
patterns  visible  in  Band  6  (0.7  to  0.8  micrometers)  or  3and  7  (0.3  to 
1.1  micrometers)  are  almost  certainly  due  to  surface  conditions  such 
as  sunglint,  foam  on  the  water,  or  turbidity  at  the  surface.  Figures 
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25,  26,  27,  and  28  show  similar  image  patterns  Ip.  all  four  spectral 
bands.  These  computer  generated  1 1  grey- map s"  were  produced  from  the 
Landsat  computer  compatible  tape  (CCT)  for  scene  (1724-14472  of 
17  July  1974). 

When  image  patterns  are  determined  to  be  the  result  of  bottom 
reflection,  it  is  usually  desirable  to  determine  the  depth  of  water 
at  that  location.  Methods  of  water  depth  determination  from  Landsat 
data  have  been  under  investigation  for  nearly  ten  years.  Work  at  the 
Environmental  Research  Institute  of  Michigan  (ERIM)  for  the  National 
Aeronautics  and  Space  Administration  (NASA) ,  the  Defense  Mapping  Agency 
(DMA),  and  the  Office  of  Naval  Research  (ONR) ,  is  encouraging.  Remote 
mapping  of  water  depth  from  multispectral  scanner  data  was  first  demon¬ 
strated  with  aircraft  data  collected  over  Lake  Michigan  (Brown,  et  al., 
1971).  The  model  developed  during  that  study  was  successfully  tested 
during  the  NASA/Cousteau  Ocean  Bathymetry  Experiment  when  depths  de¬ 
termined  from  Landsat  MSS  data  agreed  within  10  percent  of  charted 
depth  in  waters  to  42  meters  deep  (Polcyn,  1976).  Color-coded  maps  of 
the  Little  and  Great  Bahama  Banks  were  generated  with  high  gain  Landsat 
multispectral  scanner  data  showing  variations  in  water  depth  to  demon¬ 
strate  the  feasibility  of  satellite  bathymetric  mapping  (Polcyn,  Lyzenga, 
and  Tanis, .  1977).  Further  progress  was  reported  by  Lyzenga  and  Polcyn 
in  1979,  and  Denoyer  reported  successful  use  of  the  depth  determination 
algorithm  by  the  U.S.  Geological  Survey  in  analysis  of  the  Palau  Islands 


(1979) . 
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Procedures 

Multispectral  analysis  techniques  were  used  to  determine  if 
patterns  seen  over  the  shoals  were  related  to  depth.  This  was 
accomplished  by  using  a  water  reflectance  model  which  accounts  for 
reflection  of  light  from  the  ocean  bottom  developed  at  the  Environ¬ 
mental  Research  Institute  of  Michigan  (ERIM)  (Polcyn,  1976).  Digital 
numbers  calculated  from  the  model  were  compared  with  actual  digital 
numbers  produced  by  the  Landsat  MSS  to  test  whether  patterns  seen  in 
the  image  were  due  to  bottom  reflectance  or  to  surficial  phenomena. 

The  water- reflectance  model  used  has  been  tested  and  proved  useful  for 
measuring  water  depths  within  10  percent  accuracy  to  depths  of  20  meters 
in  clear  ocean  water  (Doak,  et  al.,  1979).  The  environment  of  Cape  Cod 
does  not  meet  all  the  assumptions  of  the  model,  however,  and  results 
of  the  model  are  of  limited  accuracy.  Despite  this,  use  of  this  model 
was  expected  to  aid  in  establishing  whether  or  not  features  seen  in  the 
images  can  be  correlated  with  reflected  light  from  the  sea  bottom. 


Water  Reflectance  Model  and  Estimation  of  Parameters 

The  energy  received  by  the  Landsat  MSS  sensors  in  band  (i)  over 
clear  shallow  water  can  be  computed  from: 


or,  (VrVsi) 
where : 


vi  =  vsi  +  Ksi 


V  •  e’2kiZ 
01 


t1t2 

*n^ 


Eoi  TiRbi  e'2kiZ 


=  Landsat  MSS  signal  value  in  volts 

*V  =*  Signal  over  deep  water 
si 
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Ksi  =  Landsat-1  low  gain  sensitivity  constant  (Taranik,  1978) 

*  (127/2.48)  centimeters^  steradian/milliwatt  for  MSS  Band  4 

*  (127/2.00)  centimeters”  steradian/milliwatt  for  MSS  Band  5  ■ 

*  (127/1.76)  centimeters^  steradian/milliwatt  for  MSS  3and  6 

=  (63/4.60)  centimeters”  steradian/milliwatt  for  MSS  Band  7 

T1  “  t2  “  0.98  *  Air- water  boundary  transmittance  j 

n  =  1.33  *  index  of  refraction  of  water 
*Eq^  *  Total  irradiance  in  the  spectral  band  (i)  at  the  , 

water  surface 

*  Atmospheric  radiance  transmittance  *0.8 

*  Bottom  reflectance 

*  Water  attenuation  coefficient 

i 

2  *  Depth  of  ocean  bottom  j 

j 

*Voi  *  Signal  as  depth  goes  to  0,  after  subtraction  of  deep  j 

water  signal 

The  five  variables  indicated  by  *  are  those  for  which  data  needed  to  be 
derived  or  estimated  from  other  information  to  use  the  model  without 
field  observations  in  the  study  area. 

Solar  Irradiance  at  Sea  Level 

Initially,  the  value  for  E0^  was  calculated  for  all  Landsat  MSS 
bands  from  tables  of  solar  irradiance  at  sea  level  for  an  air  mass  of 
two  (Wolfe  and  Zissis,  1978).  The  Turner  atmospheric  model  (Turner, 

1974)  was  used  to  improve  results  of  initial  calculations  from  tabular 
data  (Wolfe  and  Zissis,  1978)  by  more  accurately  modeling  the  shorter 
wavelength  content  of  the  solar  spectrum  for  a  sun  elevation  of  57°, 

! 

i 


V 
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23  kilometers  visibility,  and  an  air  mass  of  one  at  sea  level.  The 
computer  printout  for  the  solar  irradiance  is  shown  in  Figure  29. 
Improved  values  used  were: 


04  * 

9.6700 

2 

milliwatts /centime ter 

steradian 

05  * 

9.0225 

milliwatts /centime ter^ 

steradian 

Sun  elevation  57° 

06  " 

7.0125 

milliwatts/centimeter^ 

steradian 

Visibility  23 
ki  lometers 

07  “ 

16.06125 

milliwatts /centime ter^ 

steradian 

Air  Mass  of  one 

These  values  of  irradiance  (incident  to  the  normal  plane)  were 
multiplied  by  cosine  33°,  to  correct  for  the  angle  of  incidence  of 
direct  solar  irradiance  at  the  plane  of  the  water  surface. 

Deep  Water  Signal 

In  order  to  use  the  water- reflectance  model  as  written  (V^-Vs^)  * 
VQi  e'2kiZ,  one  must  first  determine  the  deep  water  signal  (Vgi)  in  all 
bands.  The  deep  water  signal  (Vs±)  (the  lowest  MSS  spectral  signal) 
represents  the  signal  when  direct  solar  irradiance  is  completely  attenu¬ 
ated. 

Eight  likely  areas  within  the  scene  were  evaluated.  The  lowest 
signals  in  all  bands  were  found  in  the  area  of  "dark  water"  east  of 
Monomoy  Island  (Figures  1  and  30).  Water  depth  at  this  location  is 
charted  as  32.4  meters  (107  feet).  Computer  listings,  showing  digital 
numbers  (signal  strength)  in  all  four  MSS  bands,  were  generated  with 
the  11  LINE  software  program  developed  at  ERIM.  This  program  uses 
Landsat  Computer  Compatible  tape  (CCT)  data  as  input  for  manipulating 
and  processing  the  data  into  various  statistical  and  visual  hard  copy 
formats.  Histograms,  statistics,  direct  output  of  wave  band  digital 
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numbers  (signal  values)  and  line  printer  "grey-maps"  can  be  generated. 
Selected  within  "dark-water"  areas,  6x6  pixel  subareas  were  averaged 
to  obtain  the  lowest  deep  water  signal.  The  deep  water  signals,  in 
digital  numbers,  were:  MSS  3and  4,  18.63888;  MSS  3and  5,  8.47222;  MSS 
Band  6,  4.1666;  MSS  Band  7,  0.38888.  The  apparent  radiance  was  calcu¬ 
lated  by  dividing  the  digital  number  by  the  sensitivity  coefficient  for 
each  MSS  Band.  Apparent  reflectance  was  calculated  using  the  formula 

R  =  11  L  ,  where:  R  =  apparent  reflectance 
Eoi 

ir  *  factor  to  convert  cone  of  radiance  to  point 
source 

L  =  apparent  radiance 

EQi  =  irradiance  of  direct  sunlight  upon  the  water 
surface . 

These  values  were  plotted  on  semilog  graph  paper,  together  with  plots 
of  other  reflectance  values  from  other  areas  within  the  scene  (Figure 
31).  The  deep  water  signal  is  the  lowest  obtained  over  all  the  water 
areas.  It  is  very  significant  that  the  deep  water  signal  was  found 
in  an  area  shown  to  be  calm,  from  direct  observation  from  the  Nantucket 
Lightship. 


Sand  and  Foam  Reflectivity 

The  solar  irradiance  and  deep  water  signal  parameters  for  the 
water  reflectance  model  were  easily  derived  from  data  on  hand.  To 
estimate  the  reflectance  of  the  sea  bottom  for  a  geographic  area  800 


miles  away  was  more  of  a  problem.  Sand  can  be  highly  variable  in 


reflectance,  and  studies  of  sea  bottom  reflectivity  are  uncommon. 
Aircraft  imagery  and  photography  over  the  study  area,  from  which 
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values  of  reflectivity  could  be  derived,  are  also  scarce.  A  study 
of  upwelling  radiance  in  the  area  for  NCAA's  National  Environmental 
Satellite  Service  (NESS)  had  been  recently  initiated  (Wright,  1979'), 
but  did  not  appear  useful  to  the  specific  problem  of  this  thesis. 

John  Portney  volunteered  to  send  the  author  a  sample  of  beach  sand 
from  Cape  Cod  National  Seashore  (1979),  from  which  the  spectral  re¬ 
flectivity  of  the  sand  could  be  measured  in  the  laboratory.  While 
one  sample  may  not  be  fully  representative  of  an  entire  pixel,  and  sand 
from  mainland  Cape  Cod  may  no't  be  representative  of  Nantucket  Shoals, 
measurements  from  this  sample  were  believed  better  than  no  measure¬ 
ment  at  all.  The  Cape  Cod  sand  sample  came  from  a  high  berm  formed 
by  Hurricane  Camile,  and  represents  the  larger  grain- size  sand 
normally  found  offshore  in  high  energy  wave  environments.  It  was 
acquired  on  30  September  1979,  on  the  east  coast  of  Cape  Cod,  near 
South  Wellfleet.  It  is  not  an  unreasonable  assumption,  considering 
the  geology  of  the  area,  that  quartz-rich  sand  eroded  from  Cape  Cod 
would  be  carried  by  littoral  drift  currents  and  deposited  upon  the 
shoals.  Mineral  composition  of  the  sample  was  probably  similar  to 
that  of  sediments  near  Nantucket.  Indeed,  Nantucket  Island,  itself, 
represents  the  glacial  moraine  sediments. 

The  effects  of  sea  state,  breakers  (cresting  waves)  and  resulting 
sea  form  should  be  a  consideration  when  evaluating  the  upwelling  re¬ 
flectance  to  the  sensor;  therefore,  to  set  a  benchmark  estimate  of 
what  sea  foam  may  look  like  radiometrically ,  a  sample  of  Gillette 
Foamy  Shaving  Cream  was  used  as  a  test  material  to  simulate  sea  foam. 
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Spectral  reflectivity  measurements  of  the  sand  and  foam  samples 
were  made  with  a  Beckman  DK-2A  Spectroref lectometer  at  ERIM  on 
19  October  1979.  Samples  of  dry  sand  and  foam  were  each  placed  in 
tinfoil- lined  film  container  cups.  Glad-Wrap  clear,  polyethylene 
film  was  used  to  cover  the  sand  sample.  3arium  sulfate  (BaSO^)  was 
used  as  the  reflectance  standard.  Instrument  output  was  a  series 
of  analog  plots  of  total  hemispherical  reflectance  over  wavelength  for 
the  spectral  band  from  (  X  ^)  to  (  x  2)  (Figure  32). 

Reflectivity  of  the  foam  and  the  film  was  read  directly  from 

the  analog  plots.  However,  because  the  reflectivity  of  the  film 

varies  from  one  wavelength  to  another  and  because  a  portion  of  the 

radiant  energy  is  reflected  a  second  time  at  the  film  plane  before 

exiting  to  be  recorded,  the  reflectivity  of  the  sand  could  not  be 

read  directly  from  the  analog  plots.  Reflectivity  of  the  combination 

of  sand  and  film  components  is  represented  by  the  following  formula: 

0  (combination  of)  =  P  (film)  +  V^fil^f^nd) 
sand  and  film 

which  was  solved  for  the  reflectivity  of  sand  at  each  wavelength 
where:  P  =  reflectivity 

T  =  Transmittance 
P  (film)  *  1  -T(film). 

Reflectivity  for  each  sample  (and  transmissivity  of  the  film) 
was  determined  at  the  midpoints  of  the  Landsat  bands  of  operation  with 
the  following  results: 
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Band 

Spectral  Midpoint 

o  Sand 

P  Foam 

t  film 

4 

5.5  micrometers 

0.28 

0.96 

0.905 

5 

6.5  micrometers 

0.36 

0.96 

0.908 

6 

7 .5  micrometers 

0.40 

0.967 

0.916 

7 

9.5  micrometers 

0.50 

0.92 

0.915 

Reflectivity  of 

the  shaving  cream  was 

nearly 

as  high  as 

the  barium 

sulfate  calibration  plate.  The  effect  of  sea  foam  covering  a  portion 
of  a  water  area  viewed  by  the  sensor  can  be  estimated  from: 

P  =  (1-F)  0  water  +  (F)  o  (foam) 

where  F  is  the  fraction  of  the  resolution  element  covered  by  foam.  For 
example,  using  data  from  MSS  Band  7,  where  the  reflectivity  of  water 
equals  0.005,  calculation  shows  that  if  foam  covered  5  percent  of  the 
water  area,  the  resulting  signal  would  be  increased  by  a  factor  of  ten, 
assuming  transmission  through  the  foam  is  negligible  and  that  real  sea 
foam  is  as  reflective  as  the  shaving  cream. 


Water  Attenuation  Coefficient 

To  estimate,  the  water  attenuation  coefficient,  the  last  required 
parameter  of  the  water  reflectance  model,  a  calibration  site  was 
selected  on  the  gently  sloping  Billingsgate  Shoal,  located  on  the  east 
side  of  Cape  Cod  Bay  where  it  is  somewhat  better  protected  from  waves 
by  the  shorter  fetch  of  Cape  Cod  Bay  (Figures  1,  4,  30,  and  33).  The 
shoal  has  a  long  and  distinctive  triangular  shape  with  a  gentle  slope 
towards  the  deeper  water  of  Cape  Cod  Bay.  Charted  depths  can  be  easily 
correlated  with  line  printer  "grey-maps"  to  obtain  line  and  point  co¬ 
ordinates  of  digital  numbers  (signal  counts)  for  the  Billingsgate  Shoal 
area. 


A  least  squares  regression  was  carried  out  to  find  the  best 
linear  fit  between  the  line  and  point  coordinates  of  digital  numbers 
(signal  counts)  and  the  geographic  coordinates  from  the  bathymetric 
charts  (Mercator  projection)  for  selected,  distinctive  control  points, 
and  yielded  the  following  equations: 

Line  =  2132.00-22.4130  (Latitude)  +  3.3988  (Longitude) 

Point  *  774.51-  8.3388  (Latitude)  -  23.0160  (Longitude) 

Files  of  positioned  depths  without  tidal  correction  were  then 
regressed  against  files  of  the  logs  of  digital  numbers  (signal  counts) 
in  MSS  3ands  4  and  5.  Tidal  information,  although  requested,  had  not 
been  received.  The  slope  of  the  regression  line  in  each  MSS  band 
represents  the  apparent  water  attenuation  coefficient  for  that  band 


of  operation.  Scatter  plots  of  the  two  bands  are  shown  in  Figures 
34  and  35.  The  values  for  VQ  and  k  were  calculated  from: 


Vi  '  Vsi  "  Voi  e 


or: 

In  (Vi-Vai)  =»  In  VQi  -  2kZ . 

The  computed  values  for  VQi  and  k  are  then: 

MSS  Band  4:  VQ  =  el. 7887  -  5.98;  k  =  0.01528  (foot*1)  = 
0.05015  (meter-1) 

MSS  Band  5:  VQ  ■  el. 6455  »  5.18;  k  =  0.2215  (foot*1)  = 
0.07267  (meter*1) 


Comparing  these  values  for  MSS  Band  4  and  MSS  Band  5  attenuation 
to  values  obtained  by  Jerlov  (1976),  reveals  that  the  MSS  3and  4  value 
of  0.0515  (meter*1)  is  about  the  same  as  the  attenuation  coefficients 


reported  for  areas  of  deep  ocean.  A  ratio  of  MSS  3and  5/MSS  Band  4 
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should  yield  a  value  of  approximately  0.25  and  values  for  MSS  Band 

4  should  range  between  0.1  -  0.2  (meters-^*)  (Lyzenga,  1980).  MSS  Band 

5  values  should  range  between  0.35  -  0.45  (meters'^*).  The  value  for 
MSS  Band  5  seemed  far  too  low  to  be  considered  reasonable,  because  0.07 
(meters'^)  is  lower  than  the  absorption  coefficient  for  water  at  that 
wave  length.  Because  of  this,  the  trwo-band  calculation  of  water  depth 
was  not  used,  but  the  single-band  calculation  seemed  to  hold  some 
promise  for  successfully  relating  the  observed  phenomenon  to  depth. 

The  maximum  depth  that  can  be  mapped  using  the  MSS  Band  4  attenuation 
coefficient  would  be  equivalent  to  one  digital  number  above  the  deep 
water  signal.  That  is*. 

maximum  depth  Z  -  depth  corresponding  to  V-Vs  *  1  count 

or,  In  (V- Vg )  *  0.  In  this  example  it  would  be 

In  (V-V  )  -  (0.7887  -  0.030564Z) 
s 

Z  -  52.624  -  28.652  In  (V-Vs) 

Z  maximum  =*  16  meters. 

Experimental  Displays  of  Tests 

A.  symbolic  display  of  depths,  by  single  band  calculation,  was 
generated  on  the  computer  for  the  calibration  site  (Figure  36).  Manual 
calculations  using  averaged  pixel  values  showed  depth  values  to  be 
within  reason  for  the  chart,  last  surveyed  in  1957. 

The  values  of  calculated  depths  seemed  to  be  within  reasonable 
(that  is  to  say  tolerable)  limits  for  the  Billingsgate  Shoal  cali¬ 
bration  site.  The  digital  display  indicated  a  depth  of  roughly  2.5 
meters  +  2.5  meters  compared  to  a  3.4  meters  plotted  depth  at  the 


computer  display  coordinates  line  934,  point  215.  The  actual  3x3 
pixel  values  compared  to  calculated  theoretical  signal  values  are  as 
follows : 


MSS  Band 

4 

5 

6 
7 


Actual  Signal  Calculated  Theoretical  Signal 


27.8889 

28.7163 

14.0000 

10.3689 

7.8889 

4.2003 

1.7778 

0.3889 

The  actual  reflectivity  plotted  with  the  calculated  theoretical  re¬ 
flectivity  is  shown  in  Figure  7 . 

A  display  of  the  Pollock  Rip  subarea  is  shown  in  Figure  38..  This 
display  may  be  compared  to  the  MSS  Band  4  "grey- map"  in  Figure  39,  and 
to  the  portion  of  the  chart  of  the  same  area  shown  in  Figure  40. 

A  depth  map  generated  for  the  Old  Man  Shoal  area  off  Nantucket 
(Figure  41),  proved  less  successful.  While  it  appears  to  show  the 
shoal  patterns,  depth  ranges  do  not  adequately  agree  with  the  charted 
depths  (Figure  42).  A  graphic  plot  of  Old  Man  Shoal  showing  the 
bottom  profile  scaled  from  the  charted  bathymetric  contours  is  com¬ 
pared  with  the  actual  MSS  Band  4  reflectivity  along  display  line 
1799  in  Figure  43.  The  geographic  position  of  this  line  is  shown 
in  Figure  42.  Actual  reflectivity  was  calculated  using  a  "moving 
window"  averaging  3x3  pixels  along  the  display  line,  then  plotting 
the  averaged  reflectance  corresponding  to  the  pixel  position  of  the 
depth  contour.  Figure  44  shows  the  plot  of  the  actual  signal,  the 
averaged  signal,  and  the  computed  theoretical  signal  to  be  expected 


for  pixel  721  in  line  1799.  The  y-Vg  reflectivity  is  plotted  in 
Figure  45  and  indicates  higher  reflectivity  in  MSS  Band  4  and  MSS 
Band  5. 

These  two  graphs  indicate  that  the  actual  signal  received  by 
the  sensor  over  the  shoal  is  not  the  same  as  the  signal  expected  from 
theoretically  calculated  values  related  to  reflection  of  light  from 
the  ocean  bottom.  The  same  calculations  were  performed  many  times 
for  other  areas  in  the  Nantucket  Shoals,  and  the  pattern  of  the 
graphic  plots  was  similar. 

A  typical  result  of  the  plot  of  actual  V-Vg  reflectivity  is 
shown  for  an  area,  (pixel  coordinates  line  1789,  point  708),  which 
the  author  called  Shoals  "A,”  where  the  water  depth  was  9  meters 
(Figure  46).  The  reflectivities  of  the  actual  signals  for  MSS  Band 
6  and  MSS  Band  7  are  much  higher  than  the  calculated  theoretical 
values  for  MSS  Band  6  and  MSS  Band  7.  Calculations  indicate  that 
the  water  attenuation  coefficients  for  these  bands  are  so  large  that 
there  would  be  negligible  penetration,  thereby  producing  a  signal 
of  zero.  V-Vg  reflectivity  plots  consistently  show  higher  reflect¬ 
ance  for  MSS  Band  4  than  the  other  bands.  Figures  47  and  48  are 
the  linear  and  log  scale  plots  of  V-Vg  reflectivity.  When  computed 
theoretical  values  of  V  were  subtracted  from  the  actual  signal  values 
in  each  band,  and  reflectivity  was  plotted,  the  plotted  line  is 
flatter  and  higher  than  zero.  This  suggests  a  contributing  factor 
independent  of  wavelength.  The  proportions  of  the  direct  solar 
irradiance  spectrum  (E  ^)  are  the  same  as  the  proportions  of  the  in¬ 
creased  radiance  (L^)  of  the  V  -  theoretical  V  shown  in  Table  5, 


therefore,  surface  specular  reflection  is  the  most  probable  contribu 
ting  agent. 


Another  peculiarity  of  the  data  set  for  this  image  is  that  MSS 
Band  7  signals  are  higher  than  those  usually  encountered  (normally 
the  range  of  0  to  1  digital  numbers).  This  scene  had  occasional 
values  as  high  as  four  digital  counts.  Landsat  data  gathered  by 
ERIM  for  the  Bahamas  test  range  did  not  exhibit  such  high  values 
CLyzenga,  1980).  Figures  49  and  50  show  the  MSS  Band  7  "grey-map" 
for  Monomoy  and  Nantucket  Islands  and  environs,  and  show  the  distri¬ 
bution  of  digital  numbers  higher  than  one  for  water  areas.  The 
higher  counts  cluster  in  areas  of  shoals  and  rips. 

Effects  of  Sunglint 

The  V-Vg  reflectivity  for  the  "grey -water"  areas  in  the  image 
was  calculated  to  see  what  the  reflectance  spectrum  looked  like. 
Sunglint  was  suspected  as  the  cause  of  the  "grey -water"  areas.  The 
familiar  pattern  of  the  spectral  signature  of  the  shoals  reappeared 
(Figure  51). 

The  similar  pattern  of  the  V-Vg  reflectivity  graphs  was  enough 
cause  to  suspect  that  sunglint  was  the  primary  factor  in  the  "grey- 
water"  reflectivity,  but  further  evidence  was  necessary.  Analysis 
of  the  wave  slope  required  to  specularly  reflect  incident  solar 
irradiance  into  the  view  of  the  sensor  also  supports  the  hypothesis. 
Considering  that  the  MSS  angle  of  view  is  5°  and  that  the  solar  ele¬ 
vation  at  the  time  of  imaging  was  57° ,  an  estimation  of  the  wave 
slope  necessary  to  reflect  direct  sunlight  in  the  direction  of  the 
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sensor  is  shown  in  Figure  52.  Also,  Che  satellite's  orbital  path 
is  approximately  10°  east  of  north,  and  at  that  time  of  year  the  sun 
in  the  east  (azimuth  of  120°)  is  more  in  the  direction  (within  10°) 
of  the  line-scan  direction  of  the  MSS  scanner  (Figure  53).  The 
probability  of  sunglint  will  increase  as  the  direction  of  scan 
approaches  the  direction  of  the  sun. 

A  much  more  rigorous  solution  to  determine  the  wave  slope 
necessary  to  reflect  solar  irradiation  "down  the  throat"  of  the 
sensor  uses  unit  vectors.  Figure  54  illustrates  this  approach  to 
the  problem,  where: 

S  *  Solar  direction  unit  vector  =  (sin  33°  sin  120°; 
sin  33°  cos  120°;  cos  33°) 

V  =  View  direction  unit  vector  =  (sin  5°  sin  100°; 
sin  5°  sin  100°;  cos  5°) 

N  =  Wave  normal  \  unit  vector  =  (S  +  V) 

(Is  +  V]) 

The  result  is  a  vector  that  has  the  coordinates: 

X(0.278) ,  Y(-0. 144) ,  and  Z(0.917) 

X 

The  orientation  of  the  wave  is  the  arc  tangent  of  —  =  1.939  =  62.7° 

The  slope  of  the  wave  is  the  complement  of  the  arc  cosine  of 

X  =  0.304  =  (90°  -  72.30°)  =  17.69° 

Of  the  three  wind  directions  recorded  for  the  time  of  imaging 
(Table  2),  both  Chatham  Coast  Guard  Station  and  Nantucket  Lightship 
recorded  a  wind  direction  from  the  north- northeas t ,  and  Nantucket 
Lightship  recorded  a  swell  direction  of  330°.  The  X  axis  for  the  compu¬ 
tation  is  east,  so  90°  should  be  added  to  62.7°  to  equal  152. 7°,  which 
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is  nearly  perpendicular  Co  Che  direccion  of  che  wind.  Therefore, 
because  che  orienCaCion  and  slope  of  che  wave  necessary  Co  specularly 
reflecC  irradiance  co  che  sensor  agrees  wich  che  known  condicions, 
che  chances  of  sunglinC  are  nearly  cercain.  The  conclusive  evidence 
of  sunglinC,  ChaC  is  a  direcCional. change  due  Co  a  change  of  viewing 
angle  recorded  by  an  inscancaneous  sensor  (che  RBV  sensor)  is  absent. 
The  RBV  was  not  used  for  the  scenes  imaged  by  che  MSS  sensor. 

The  reflectivity  graph  was  approximately  equal  in  all  bands 
except  MSS  Band  4.  Occasionally  che  reflectance  in  MSS  Band  4  was 
as  much  as  1.6  times  higher  than  the  ocher  three  bands.  Figure  51 
is  included  to  illustrate  the  pattern  of  the  graphs.  This  phenomenon 
is  very  interesting,  because  if  specular  reflection  (sunglint)  is  the 
cause  of  che  "grey-water"  areas,  reflectivity  in  all  four  bands  must 
be  equal.  Turner's  atmospheric  model  was  used  to  recalculate  the  re¬ 
flectance  to  more  accurately  model  the  shorter  wavelength  portion  of 
the  spectrum  for  an  air  mass  of  one.  The  result  of  the  calculations 
did  not  remove  the  MSS  Band  4  "spike"  in  the  graphs. 

An  error  in  radiometric  calibration  could  cause  such  a  variation 
between  the  spectral  bands.  Studies  of  calibration  for  Landsat-1  are 
currently  in  progress  (Chance,  1980).  Errors  in  radiometric  cali¬ 
bration  have  been  noted  in  the  past  for  Landsat-2  for  the  period  of 
22  January  co  15  July  1975  by  Colvocoresses  (1977),  and  corrections 
have  been  reported  by  NASA  (Anonymous,  1977). 

However,  on  the  assumption  that  the  radiometric  calibration  for 
Landsat-1  during  the  period  in  question  is  correct,  there  are  other 
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possible  explanations  of  this  phenomenon.  White  caps  may  reflect  an 
increased  amount  of  skylight,  in  addition  to  direct  sunlight,  into 
the  field  of  view  of  the  sensor.  The  proportion  of  skylight  is 
greater  in  MSS  3and  4  and  would  add  to  the  apparent  reflectance.  If 
such  turbulence  was  present,  MSS  Band  4  would  be  more  affected  by 
the  light  scattered  from  submerged  air  bubbles  because  of  its  deeper 
water  penetation,  than  in  the  other  bands  (Lyzenga,  1980).  It  is 
reasonable  to  assume,  but  cannot  be  proven  from  the  data  available, 
that  wind  speed  may  be  higher  in  "grey- water"  areas.  A  higher  wind 
speed,  acting  in  conjunction  with  the  longer  fetch  of  the  ocean, 
could  produce  more  whitecaps  than  were  recorded  at  Chatham  (Table  2) . 
The  phenomenon  of  variation  of  reflectance  in  the  image  is,  in  itself, 
perhaps  the  best  evidence  of  higher  sea  state  conditions  when  there 
is  no  other  contrary  evidence. 

It  should  be  emphasized  that  the  water  reflectance  model  used  in 
this  study  deals  with  the  reflectance  of  light  from  the  ocean  bottom, 
and  does  not  include  the  phenomenon  of  variation  in  surface  and  near¬ 
surface  reflectance  (Polcyn,  1976).  The  signature  used  for  Figure  51 
came  from  an  area  which  was  charted  to  be  55  meters  deep. 

Another  possible  explanation  which  was  only  briefly  considered  as 
a  cause  of  the  additional  upwelling  radiance  in  MSS  Band  4  was  that 
of  marine  organisms,  such  as  phytoplankton  which  might  be  a  factor; 
however,  any  added  contribution  to  upwelling  reflectance  caused  by 
phytoplankton's  Chlorophyll  A  would  cancel  in  the  equation  used.  The 
very  small  contribution  (Figure  55)  from  phytoplankton  would  be  diffi¬ 
cult  to  detect  with  the  Landsat  MSS,  even  with  a  ship  on- sine  to 


verify  the  organism's  presence  (Hovis,  1980).  it  is  also  unlikely 
that  they  would  have  any  correlation  to  the  patterns  in  the  image. 

Turbidity  Effects 

An  area  (Figure  56)  within  Nantucket  Sound,  which  was  thought 
to  contain  turbid  water  (from  Multitemporal  Analyses),  was  evaluated 
multispectrally  to  see  how  the  plotted  signature  differed  from  other 
plots  (Figure  57).  The  slope  increased  from  MSS  Band  6  to  MSS  3and  7 
for  the  V-V3  plotted  signal  values  (compare  to  Figure  48),  The  other 
three  curves  on  Figure  57  represent  the  simulated  V-Vg  curves  compu¬ 
ted  from  the  Quasi- Single- Scattering  (QSS)  model  of  Gordon  (1973) 
using  different  water  attenuation  coefficients  for  the  Jerlov  (1976) 
water  types.  Theoretical  modeling  studies  by  Maul  and  Gordon  (1975) 
have  shown  that  the  water  reflectance  can  be  approximated  by  the 
equation: 

R  =  K  ab _ 

a  +  b 

where:  R  •  reflectivity 

K  =  geometric  factor  equal  to  0.27 
a  =  the  absorption  coefficient  of  water 

(including  dissolved  and/or  suspended  material) 
b  »  the  backscattering  coefficient  (the  fraction  of  the  total 
scattering  coefficient  which  results  in  scattering  in  the 
backward  hemisphere) 

For  these  computations  R  was  set  equal  to  the  reflectivity  of 
MSS  Band  4;  a  was  set  equal  to  Jerlov' s  coefficient  for  water  type  9 
for  coastal  water;  the  equation  was  then  solved  for  b.  The  value  for 
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fa  was  then  used  for  all  other  bands  of  operation,  assuming  that 
scattering  was  independent  of  wavelength.  The  same  procedure  was 
used  for  other  water  attenuation  coefficients  for  each  of  the  other 
water  types. 

The  graphs  on  Figure  57  indicate  that  turbidity  is  present  and 
that,  as  expected,  the  initial  constants  selected  for  use  with  the 
water-  reflectance  model  would  cause  errors  in  depth  measurement 
(Polcyn,  1976).  In  addition,  there  were  probably  surficial  effects 
which  would  also  cause  errors,  just  as  if  there  were  a  variation 
of  turbidity  density. 

The  Old  Man  Shoal  V-Vg  reflectivity  (Figures  45  and  58)  was 
examined  for  turbidity,  with  the  knowledge  that  evidence  of  sunglint 
had  been  demonstrated.  The  reflectivity  bias  at  all  four  wavelengths 
was  set  at  2.3  percent,  because  the  reflectivity  in  MSS  3and  6  and 
MSS  Band  7  was  approximately  flat  and  equal  to  2.3  percent.  The 
two  main  water- penetrating  bands  (MSS  Bands  4  and  5)  would  therefore, 
have  contributions  not  only  from  depth  but  also  from  suspected 
turbidity . 

QSS- simulated  reflectance  curves  were  generated  in  the  same  manner 
as  before,  except  that  the  bias  had  been  subtracted  from  the  MSS  Band 
4  reflectivity.  The  actual  V-Vg  curve  fell  between  the  theoretical 
curves  for  0.19  and  0.12  (k^ .meter  S  and  may  be  more  appropriate  for 
water  depth  measurements  over  this  local  area  of  the  Old  Man  Shoal, 
where  turbidity  and  sunglint  both  are  present. 
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Sea  Foam  Consideration 

Sea  foam  as  a  Lambertian  reflector  could  raise  the  signal  value 
in  the  same  manner  as  sunglint,  except  sea  foam  reflectance  is  not 
directional.  When  the  transect  across  Old  Man  Shoal  was  mapped  in 
Figure  43,  it  was  noted  that  the  highest  signal  values  did  not  fall 
exactly  on  the  peak  of  the  charted  shoal,  but  were  offset  in  the 
westward  direction  by  one  pixel.  This  could  be  explained  by  the 
action  of  waves  breaking  over  the  shoal  creating  and  spreading  sea 
foam.  Another  explanation  is  that  the  shoal  could  have  migrated  or 
changed  shape  since  it  was  last  surveyed  23  years  ago. 

Surface  Manifestations  of  Internal  Waves 

Figure  59  is  a  "grey- map"  of  the  surface  manifestation  of  in¬ 
ternal  waves.  There  are  approximately  seven  waves  across  the  display 
of  124  pixels  from  points  625-749.  There  is  also  some  apparent  re¬ 
fraction  with  other  waves  in  the  display.  The  leading  wave  is 

propogating  westward  in  the  convex  direction  of  the  arc.  The  period 
between  successive  waves  diminishes  from  approximately  1,368  meters 
at  the  leading  wave  to  approximately  625  meters  at  the  eastern  edge 
of  the  display.  The  wave  is  oriented  about  10°  east  from  north. 

A  plot  of  the  digital  numbers  in  each  multispectral  band  across 
the  leading  wave  in  line  1031,  points  625  through  654,  is  shown  in 

Figure  60.  The  wave  is  sinusoidal  in  each  band.  The  slopes  of  lines 

drawn  from  the  highest  signal  in  each  band  to  the  lowest  signal  in 
each  band  are  different  on  each  side  of  the  wave  form.  Table  6  lists 
the  approximate  slope  for  each  band  of  operation  on  the  eastern  and 
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western  sides  of  the  wave  form.  The  rate  of  signal  increase  in  each 
band  on  the  eastern  side  of  the  wave  form  is  approximately  half  that 
of  the  rate  of  signal  increase  on  the  western  side  of  the  wave  form. 

The  lowest  digital  numbers  lie  clustered  on  the  eastern  side  of  the 
plot.  The  lowest  MSS  Band  4  digital  number  is  approximately  five  digi¬ 
tal  numbers  higher  than  the  deep  water  signal  at  18.6,  so  there  is  some 
additional  radiance  other  than  diffuse  radiance  in  the  recorded  signal. 
This  observation  of  the  surface  manifestation  of  internal  waves  seems 
to  match  other  observations  in  nearly  every  respect,  even  to  the  time 
of  year,  direction  of  propogation,  and  dark  band  trailing  the  leading 
edge  of  the  wave  form.  The  shape  of  the  leading  wave  (the  actual 
physical  wave),  is  probably  different  than  the  shape  of  the  trailing 
waves,  as  suggested  by  the  pattern  of  the  wave  propagation  of  the  in¬ 
ternal  wave  itself.  This  fact  is  illustrated  in  the  work  of  Neumann 
and  Pierson  ((1966')  (Figure  38))  and  in  studies  by  LaFond  and  Cox  ((1962) 
(Figure  57))  indicating  that  the  leading  wave  has  a  higher  amplitude 
Chan. succeeding  waves.  Specular  reflection  from  the  leading  wave  sur¬ 
face  may  actually  enter  a  water  surface  nearby,  depending  on  the  shape 
of  the  wave,  and  therefore  be  absorbed.  This  could  explain  the  dark 
trailing  band  when  propogating  westward.  If  the  surface  was  slightly 
rough,  the  reflection  from  small  facets  could  explain  the  higher  signal 
received  at  the  sensor  rather  than  the  expected  lower  signal  from 
diffuse  deep  water. 

Summary  of  the  Multisnectral  Analyses 

The  water- ref lectance  model  developed  by  Polcyn  (1976)  at  ERIM, 
was  applied  to  the  low-gain  Landsat  CCT  for  the  17  July  1974  image 


(1724-14472)  of  Cape  Cod  in  order  to  aid  the  interpretation  of  hydro- 

i 

graphic  features.  The  techniques  used  to  extract  bathymetry  remotely 
from  successful  tests  in  clear  ocean  waters  were  applied  to  subareas 
in  the  scene  whose  characteristics  were  unknown.  The  purpose  was  not 
solely  to  map  depths  with  high  accuracy,  but  rather  to  determine  if 
the  signals  received  by  the  Landsat  sensor  were  from  bottom  reflectance. 

Estimates  of  the  variable  parameters  of  the  water  reflectance 
model  were  derived  from  existing  data  rather  than  through  direct 
measurement.  Transformation  of  pixel  coordinates  to  geographic  co¬ 
ordinates  was  accomplished  by  a  least  squares  regression.  Charted 
depths  and  positions  were  regressed  to  provide  an  apparent  coefficient 
of  attenuation  for  the  water  over  Billingsgate  Shoal.  Estimates  of 
the  variable  parameter  of  bottom  reflectance  and  possible  surficial 
sea  foam  were  taken  from  reflectivity  and  transmission  spectra  measured 
from  Cape  Cod  sand  and  simulated  sea  foam. 

Computer- generated  depth  maps,  using  a  single-band  algorithm  for 
depth  measurement,  and  "grey-maps"  were  created  for  the  initial  cali¬ 
bration  area  and  other  areas  within  the  image.  The  relationship  to 
depth  was  judged  by  comparison  with  charted  depths  and  comparison  with 
theoretically  computed  depths  and  theoretical  signals  calculated  from 
such  depths . 

The  signals  for  image  patterns  over  the  area  of  Old  Man  Shoal  did 
not  correlate  well  with  bottom  reflected  energy.  The  other  areas  tested 
were  more  accurate,  when  compared  to  charted  d  a-  and  to  theoretical 
signals  calculated  from  such  depth. 


Several  working  hypotheses  were  tested  to  account  for  the  cause 
of  the  errors.  Sunglint  (specular  reflectance)  was  shown  to  be  the 
primary  cause  of  a  significant  bias  in  the  signal  and  also  the  cause 
of  variations  of  reflectivity  in  deep  ocean  water  areas.  An  in¬ 
teresting  variability  of  the  reflectivity  between  the  four  MSS  spectral 
bands  of  the  MSS  sensor  was  also  observed.  This  phenomenon  was 
attributed  to  surface  turbulence  of  whitecaps,  multiple  scattering 
of  submerged  bubbles,  or  additional  reflection  of  skylight  in  the 
field  of  view  of  the  sensor. 

The  hypothesis  that  turbidity  was  present  in  the  local  water 
above  Old  Man  Shoal  was  evaluated  by  comparing  reflectances  from  the 
depth  algorithm  to  the  QSS  model  for  turbid  water  and  to  graphs  of 
calculated  theoretical  values.  Techniques  comparing  simulated  re¬ 
flectivity  curves  to  actual  calculated  reflectivities  were  used  to 
estimate  more  appropriate  water  attenuation  coefficients. 

Sunglint  and  turbidity  contribute  to  the  increased  visibility  of 
the  feature  in  question.  Digital  displays  and  observations  on  the 
identified  surface  manifestations  agree  with  previous  studies. 
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SUMMARY  AND  CONCLUSIONS 

The  Cape  Cod  study  site  represents  a  unique  marine  environment. 

It  is  subject  to  prevailing  winds  from  the  southwest  in  summer  and  to 
gales  from  the  northeast  in  winter.  The  geologically  young  (Pleistocene) 
glacial  deposits  are  subject  to  high  rates  of  erosion  and  aggradation 
from  the  action  of  winds,  waves,  and  currents.  The  physiography  of  the 
coast  and  continental  shelf  is  the  result  of  geologically  recent  innun- 
dation  from  melting  glaciers.  Offshore  sediments  represent  relict 
beaches  and  the  waters  are  relatively  clear  because  areas  of  turbid 
water  are  contained  within  estuaries,  or  the  turbid  water  is  flushed 
into  the  deep  ocean  by  the  rapid  action  of  wind,  waves,  and  strong 
currents . 

Space  imagery  can  provide  highly  detailed  information  about  the 
coastal  and  offshore  features,  some  of  which  will  never  appear  on 
standard  nautical  charts.  Bottom  areas  in  the  near  shore  marine  en¬ 
vironment  can  change  slowly  over  time  and  also  change  rapidly  in 
response  to  catastrophic  storms.  In  either  event,  the  resulting 
ocean  bottom  changes  can  represent  hazards  to  navigation.  Charts  of 
such  areas  will  require  frequent  correction  and,  if  resources  are 
available,  will  require  frequent  revision.  Detection  of  these  hazards 
to  navigation  and  their  subsequent  publication  to  the  maritime  commun¬ 
ity  is  the  subject  of  great  international  concern. 

Findings  of  this  study  support  use  of  Landsat  imagery  to  ascertain 
bottom  changes,  coastal  changes  and  man  made  changes  for  the  following 
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purposes:  Initiation  of  appropriate  marine  warnings  such  as  Notice 
to  Mariners  for  hazards  to  navigation;  and,  when  supported  by  ground 
truth  data,  chart  revision. 

Large  areal  coverage,  repetitive  coverage,  and  multispectral 
qualities  of  Landsat  data  can  provide  significant  benefits  to  the 
hydrographic  community.  Information  can  be  derived  remotely  which 
amplifies,  complements,  or  supplements  studies  done  from  aircraft 
and  ships.  Information  can  be  obtained  that  is  beyond  the  capability 
of  single  band  sensor  systems  or  from  conventional  stereo- compilation 
methodology.  Landsat  MSS  imagery  provides  additional  information  and 
furnishes  evidence  to  support  conclusions  about  observed  phenomena. 

For  example,  one  result  of  this  study  was  to  show  that  data  significant 
to  hydrographic  mapping  for  the  Cape  Cod  area  was  surprisingly  difficult 
to  obtain  even  though  the  water  around  Cape  Cod  is  an  active,  offshore 
area  for  both  economic  and  defense  purposes .  Most  space  imagery 
available  is  neither  cloud  free,  nor  without  some  sea  state.  Hydro- 
graphic  survey  data  is  out  of  date  and  oceanographic  data  is  sparse 
because  of  data  collection  hazards  posed  by  the  shoals.  With  the 
significance  of  the  Nantucket  Shoals  area  to  the  maritime  community, 
and  its  well  established  reputation  as  a  grave  yard  for  ships,  the  lack 
of  a  detailed  environmental  data  base  or  even  a  plan  to  collect  such 
data  seemed  peculiar. 

To  aid  interpretation,  relatively  new  techniques  were  applied  to 
determine  bathymetry  from  the  image  data.  Some  image  patterns  were 
identified  as  being  the  surface  manifestation  of  internal  waves. 
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Other  patterns  seen  over  the  shoals  and.  Old  Man  Shoal  in  parti¬ 
cular  (Figure  LA),  resulted  from  both  surface  phenomena  and  bottom 
reflectance.  The  Landsat  signals  received  give  erroneous  depth  measure¬ 
ments  due  to  turbidity,  and  the  signals  appear  to  be  primarily  due  to 
surface  effects.  From  all  evidence  gathered,  these  surface  effects 
include  both  sunglint  and  surface  turbulence. 

Analysis  of  "grey-water"  patterns  shows  that  their  reflectance 
spectrum  is  similar  to  those  for  shoal  areas.  All  evidence  gathered 
suggests  that  the  cause  of  the  lighter  patterns  is  sunglint  and  sur¬ 
face  turbulence;  but,  because  these  patterns  occur  in  deep  water  they 
would  not  endanger  ships. 

The  author  believes  that  this  report  demonstrates  the  need  for 
and  feasibility  of  extracting  bathymetry  data  from  Landsat  records, 
and  provides  some  guidance  to  other  users  of  these  same  data  analysis 
techniques.  The  study  of  the  sea  as  a  specular  and  diffuse  reflector 
is  a  long  term  continuing  process.  Efforts  should  be  started  to  deter¬ 
mine  the  best  way  of  estimating  water  turbidity  on  satellite  images, 
including  an  analysis  of  shadows  from  clouds  or  contrails  and  their 
utility  in  accurately  measuring  turbidity.  A  second  area  of  study 
should  be  directed  at  determining  whether  or  not  shoals  can  be  detect¬ 
ed  from  local  turbidity  densities,  as  suggested  by  the  analysis  of  Old 
Man  Shoal  in  this  paper. 


Figure  1.  Cape  Cod,  Massachusetts  and  Environs. 

(Note:  Data  for  each  image  is  provided  and  descriptors  will 
not  be  repeated  on  succeeding  caption  pages.  Data  are  as 
follows : 

Image  Identification  Number  =  1724-14472 

Sensor  and  Spectral  Band  =  Landsat- 1  band  5 

Acquisition  Date  =  17  July  1974 

Scale  =  1:1,000,000 

Wind  Speed,  Direction  and 

Reporting  Station  *  14  knots  20°  (Nantucket  Tower) 

Nantucket  Tide  Level  above 
Low  Mean  Sea  Level  (if 
Available)  =  2.3  meters 
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figure  1.  Cape  Cod,  Massachusetts  and  Environs: 
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Figure  2.  Glacial  Lcbes  ana  Deposic.cn 
Lines  (from  Scrahler,  1966). 
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Figure  5.  Color  Print  from  'Transparency  Taken 
by  Skylab  S-190B  Showing  r.hoals  of  Kantucket. 
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Figure  5-  Portion  of  Landsat-3  RBV  Image  Showing  Details  of 

Cuspate  Shoreline  in  Nantucket  Harbor  and  Subaqueous 
Hydrographic  Features. 
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Figure  5.  Portion  of  Landsat-5  REV  I sage 
Showing  Details  of  Cuspate  Snoreline  in 

Nantucket  Harbor  and  Subaqueous  Hydrographic 
Features. 


Figure  6.  btlerged  Fort  Ion  of  Ntt  lontf  J  of  Unseat  Scene  1724-1447?.  of  17  July  1974.  Bathymetry  Chart  Over 
lay  ltv<vralo  Similarity  in  Fatterna  0»«  Ihoata. 


Figure  7.  Landsat  MSS  Band  4  Image  Showing  Edge  Definition  of 
Pattern  Related  to  Direction  of  Wind. 
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Figure  8.  Landsat  MSS  Band  4  Image  Showing  Patterns  of  Sea  State. 

(Note:  Compare  water  surface  patterns  to  Figures  7-10.) 
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Figure  9.  Landsat  MSS  Band  4  Image  Showing  Few  Visible  Patterns 
in  the  Ocean  due  to  Calm  Sea  State,  High  Tide,  and 
Low  Sun  Angle.  Difference  in  Brightness  of  Land 
Surfaces  is  Attributed  to  Snow  Accumulation. 
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Figure •  Water  Maas  Boundaries  that  are  #>osaibly  L-e  to 
BiolOfeleal  Activity  are  Recorded  in  tnia  voice  alro 
Uluatr.tea  the  Effects  of  -OKl.oaetric  mSllp§>iMi£‘.m 


Figure  10.  Water  Mass  Boundaries  that  are  Possibly  due  to 


Biological  Activity  are  Recorded  in  this  Image  which 
also  Illustrates  the  Effects  of  Radiometric  "Clipping 
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Figure  11.  Coastal  Changes  of  Cape  Cod  Monitored  by  Landsat 


(From  U.S.  Geological  Survey). 
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Figure  12.  Landsat  MSS  Band  7  Image  Showing  Breach  of  Monomoy 
Island  Resulting  from  the  31izzards  of  1978, 
Deposition  from  Littoral  Drift  Currents,  and 
Variations  of  Land  Brightness  Attributed  to  Snow 
Accumulation. 
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Figure  14,  Lancia* t  HaS  Band  Intake  Showing 
•  Snail  Plune  at  the  End  of  Cane  Cod  Canal. 
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Figure  L4.  Landsat  MSS  Band  4  Image  Showing  Small  Plume  at  the 
End  of  Cape  Cod  Canal. 
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Figure  15.  Image  Showing  Sediment  Plumes  Resulting  from  Con¬ 
vergence  of  Longshore  Currents. 
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Figure  l5 .  Image  Showing  Cediaont  riuaee  Keaultir.fi 
Troa  Convergence  or  Longshore  Currents. 
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Figure  16.  Image  Showing  Shadow  of  Contrail  on  the  Water  Which 
May  Indicate  Presence  of  Turbidity. 
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Figure  17.  Image  Showing  Two  Distinct  Cloud  Types. 
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Figure  18.  Lands'*  t  MSS  linage  Showing  Major 
Patterns  in  Band  7*  (Note:  This  inuicates^ 
the  patterns  above  are  apparently  surficial' 
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Figure  18.  Landsat  MSS  Image  Showing  Major  Patterns  in  Band  7. 
(Note:  This  indicates  the  patterns  above  are 
apparently  surficial.) 
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Figure  19.  Seasat  Synthetic  Aperture  Radar  (SAR)  Image  Showing 


"Shoal  Pattern"  Visible  in  Microwave. 

(Note:  L-Band  wave  length  is  (23  cm)  and  does  not 
penetrate  water.) 

(No  image  number) 
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Figure  19.  Seasat  Syntnetic  Aperture  naaar  (oAn) 
Image  Snowing  "Snoal  ^scterr."  Visiale  ir.  'j.cro»ave. 
(i.ote:  L-iana  w_»ve  length  is  £3  ca  ana  does  not 

penetrate  water. ) 
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Figure  20 .First  of  Two  Sequential  Sky lab- 3  S-190A 
Fnctocruphs  Illustrating  the  I fractional  Cnaracterietic 
•f  Surficial  Specular  Kef lecti  or.  »r jnglir.t ; .  (Kote: 
pattern  tc  second  p;..  tocra;h* Figure  21) 
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Figure  20.  First  of  Two  Sequential  Skylab-3  S- 190A  Photographs 
Illustrating  the  Directional  Characteristic  of 
Surficial  Specular  Reflection  (Sunglint) . 

(Note:  Compare  pattern  to  second  photograph- Figure  21.) 
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Figure  21.  Second  Sequential  Skylab-3  S- 190A  Photograph 

Illustrating  the  Directional  Characteristic  of 
Surficial  Specular  Reflection  (Sunglint) . 

(Note:  The  pattern  changes  with  the  angle  of  view.) 
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Fifiure  21.  Second  Sequential  ."kv  lab-3  Sfl90A  ^hotrcrajj 
Illustrating  the  Directional  Characteristic  of  f.urljcial 
Specular  Heflecticr.  CSur^lint.)  (Kote:  Tne  patten, 
cnan&es  with  the  ar.gle  of  view.) 
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Figure  22.  First  of  Three  Sequential  Landsat  MSS  Band  j> 
Scenes  Snowing  no  Change  of  Pattern  with  Change  in  View 
Angle.  (Note:  A  charge  of  pattern  would  provide 
evidence  of  sunglint^.  see  Figures  2 3  and  2<f). 
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Figure  22.  First  of  Three  Sequential  Landsat  MSS  Band  5 

Scenes  Showing  No  Change  of  Pattern  with  Change 
in  View  Angle . 

(Note:  A  change  of  pattern  would  provide  evidence 
of  sung lint,  see  Figures  23  and  24.) 
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Figure  23-  Second  of  Three  Sequential  Landsat  MSS  Band  5  Scenes 

Showing  Mo  Change  of  Pattern  with  Change  in  View  Angle. 
(Note:  A  change  of  pattern  would  provide  evidence  of 
sunglint,  see  Figures  22  and  24.) 
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Figure  2J.  Second  of  Three  Sequential  Landsat  MS£  Band  p 
Scenes  Showing  no  Change  of  Pattern  with  Caange  in  View  Angle* 
(tote:  A  change  of  pattern  would  provide  evidence  of  sunglir.t 
see  Figures  and  Sh), 
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7i|ttrc  24.  Third  of  Three  Sequential  Lendeat  NS£  Band  ? 
Scenes  Showing  no  Change  of  Pattern  with  Change  in  View 
Angle.  (Note:  A  change  of  pattern  would  provide  eviaer.ce 
of  Aung lint,  see  Figures  22  and  23). 
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Figure  24.  Third  of  Three  Sequential  Landsat  MSS  Band  5  Scenes 

Showing  No  Change  of  Pattern  with  Change  in  View  Angle. 
(Note:  A  change  of  pattern  would  provide  evidence  of 
sung  lint,  see  Figures  22  and  23.) 
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Figure  27. "Grey- .'lap"  M3S6  Old  Man  Shoal. 
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Figure 28.  "Grey- Map"  Mo37  via  .'Ian  or.oal. 
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TABLE  1.  (coat 


TABLE  2.  SELECTED  WEATHEK  DATA  FOB  DATES  0?  IMAGE  ACQUISITION 
(see  Table  1) 

SHELL  WIND  WAVES 

Local  Visibility  Period  Height  Period  Height 

Date  I -a  l  N  long  U  Tlae  Direction  Speed  N.  Miles  in  Seconds  In  Matera  Direction  in  Seconds  In  Metera 


TABLE  3.  TIDAL  SUMMARY  FOR  NANTUCKET  ON  DATES  OF 
IMAGERY  ACQUISITION  (See  Table  1) 


Date 

LST** 

Height  in  Feet 
Above  Mean  Low  Tide 

13  Dec  73 

0953 

Flood 

4.9 

(1.5  m) 

17  Jul  74 

0947 

Ebb 

7.6 

(2.3  m) 

22  Apr  75 

0944 

Ebb 

7.1 

(2.2  m) 

28  May  75 

0944 

Flood 

4.4 

(1.3  m) 

6  Nov  76 

0942 

Flood 

6.3 

(1.9  m) 

13  May  76 

0912 

Ebb 

6.4 

(1.9  in) 

11  Apr  77 

0926 

Flood 

2.5 

(0.8  m) 

9  Mar  78 

1104 

Flood 

7.8 

(2.4  m) 

15  Apr  78 

0943 

Ebb 

6.3 

(1.9  m) 

19  Aug  78 

0944 

Flood 

9.3 

(2.8  m) 

*  Time  of  tides  calculated  from  Observatory  Station  Data,  Boston, 
Massachusetts 


**  Local  Standard  Time 


TABLE  4.  TIDAL  CURRENT  DATA  AT  POLLOCK  RIP  CHANNEL  ON  DATES  OE  IMAGE  ACQUISITION 

(see  Table  1) 


TABLE  5.  PERCENTAGE  OF  INCREASED  RADIANCE  COMPARED  TO 
SOLAR  IRRAD  LANCE 


Landsat 

MSS 

Band 

V-Vg(L) 

Radiance 
nr  /cm^sr 

2  CV-Va) 
Radiance 

CL) 

CE0) 

Irradlance 

mw/cm^ar 

X  CE0) 
Irradlance 

4 

.0848242 

22.430 

8.1099 

22.777 

5 

.0811401 

21.455 

7.5669 

21.252 

6 

.0710987 

18.800 

6.4588 

18.140 

7 

.1411167 

37.315 

13.470 

37.831 

Total : 

.3781797 

100 

35.6056 

100 
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TABLE  6. 

SLOPE  OF  DIGITAL  NUMBER  TO 
FOR  SURFACE  MANIFESTATIONS 

LINEAR  PIXEL  POSITION 
OF  INTERNAL  WAVE 

Landsat 

MSS 

Westward 

Eastward 

Band 

Slope 

Slope 

4 

35° 

20° 

5 

31° 

16° 

6 

34° 

16° 

7 

20° 

10° 
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